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We study a class of quantum measurement models. A microscopic object is entangled with a macro- 
scopic pointer such that each eigenvalue of the measured object observable is tied up with a specific pointer 
deflection. Different pointer positions mutually decohere under the influence of a bath. Object-pointer 
entanglement and decoherence of distinct pointer readouts proceed simultaneously. Mixtures of macroscopi- 



cally distinct object-pointer states may then arise without intervening macroscopic superpositions. Initially, 
object and apparatus are statistically independent while the latter has pointer and bath correlated accord- 
ing to a metastable local thermal equilibrium. We obtain explicit results for the object-pointer dynamics 
> 

depend on details of the pointer-bath coupling if it is smaller than the bath correlation time, whereas in the 



with temporal coherence decay in general neither exponential nor Gaussian. The decoherence time does not 



opposite Markov regime the decay depends strongly on whether that coupling is Ohmic or super-Ohmic. 
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I. INTRODUCTION 

The interpretation and theoretical description of measurements on quantum systems have been 
under intense debate since the birth of quantum mechanics |l|, |2|. In the last three decades the 
major role played by environment-induced decoherence in a measurement process has been fully 
acknowledged thanks to the works of Zeh, Zurek, and others (see [3j, |4j, |5|, Ig] and references therein). 
A renewal of interest for quantum detection and decoherence came in the last decade with new 
developments in quantum information. It is desirable to better understand the relation between 
quantum and classical information and how one can convert one into another. Moreover, good 
control over all sources of decoherence is required for quantum information processing. On the 
experimental side, measurements can be used either to extract information on quantum states or 
to monitor q uantu,„ SyS ,e mS Quantum trajectories fl, q ua n ,„ m Ze»„ effect QQ). Experimental 
data are now available for the decoherence time in microwave cavities |10l, in trapped ions IllH . in 

n n 

id state devices like quantum dots [TJ] and superconducting tunnel junction nanocircuits [13, 



so 

14 1 . for fullerene molecules decohered by collisions with a background gas lffl, and for beams of 

electrons decohered by Coulomb interaction with a semiconducting plate 16|. These and other 



experiments call for studies of concrete models for quantum measurements. Various models have 



been investigated so far (see, e.g., [Tjj and the interesting statistical physics models of Refs. [18J, 



3) 



but a satisfactory treatment of decoherence resulting from the many-body interactions in the 
measurement apparatus is still lacking. 

A measurement on a quantum system consists in letting this system (called "object" in the 
following) interact with a measurement apparatus, in such a way that some information about the 
state of the object is transferred to the apparatus. As already recognized by Bohr, even though the 
composite (object and apparatus) system has to be described by quantum theory, some part of the 
apparatus (called the "pointer" in the following) must be capable of classical behavior. The inter- 
action builds up a one-to-one correspondence between the eigenvalues s of the measured observable 
S (supposed here to have a discrete spectrum) and macroscopically distinguishable pointer states 
(characterized e.g. by sharply defined pointer positions separated by macroscopic distances). In 
addition to this object-apparatus coupling, the measurement must involve some "super selection 
rules" destroying the coherences between the pointer states [3j, |4J, la, |2(|. Most previous discus- 
sions in the literature consider these two processes separately: A first step ( "premeasurement" ) 
exclusively treats the unitary evolution entangling object and pointer. For an object initially un- 
corrected to the apparatus and in a linear superposition \ip$) = ^2 s c s \s) of eigenstates of S, this 



entanglement produces a superposition of macroscopically distinguishable object-pointer states 
l^cnt) = ^2 S c s \s) ® \ip-p), where \ipp) is the pointer state tied up with the eigenvalue s (for simplic- 
ity we provisionally assume that this state is pure). The "Schrodinger cat" state l^cnt) is taken as 
the initial state for a second dynamical process, decoherence, which leads to superselection rules; 
there, the quantum correlations between object and apparatus are transformed into classical cor- 
relations, as the superposition of object-pointer states is degraded to a statistical mixture of the 
same states according to l^cntX^entl — ► Yls l c «| 2 \ s )( s \ ® IV'-p)(V ; pl- F° r such a sequential treatment 
to make physical sense, the duration of the entanglement process must be short compared with the 
decoherence time idee associated with the latter transformation. However, it is known that idee is 
extremely short for macroscopic superpositions. The present paper is devoted to the more realistic 
situation where entanglement and decoherence proceed simultaneously. If the characteristic time 
for entanglement is larger than £<jec> macroscopic superpositions decohere to mixtures faster than 
entanglement can create them. The measurement process then yields the final mixture without 
involving a Schrodinger cat state at any previous moment, the object-pointer initial product state 
being directly transformed as 

mm ® i^x^i — Y. \c s \ 2 \s)( s \ ® i^x^i . (i) 

s 

Here \ips) = ^2 s c s \s), \ip^>), and |^-p) refer to the object initial state, the pointer initial state, and 
the pointer state tied up with s, respectively. 

Our model is a three-partite model and consists of the quantum object to be measured, a single 
"pointer" degree of freedom of the apparatus singled out by its strong coupling to the object and 
by affording a macroscopic range of "read-out" values, and a "bath" comprising all other degrees 
of freedom of the apparatus. A pointer-bath coupling is responsible for decoherence. Overcoming 
limitations of many previous approaches, we (i) allow for object-pointer entanglement and deco- 
herence of distinct pointer readouts to proceed simultaneously, (ii) cope with initial correlations 
between pointer and bath by considering them initially in a metastable local thermal equilibrium, 
and (iii) go beyond the Markovian treatment of decoherence. The physical relevance of point (i) 
has been discussed above. This simultaneity of entanglement and decoherence and the possibility 
of having decoherence much faster than entanglement have been considered in 18U19J]. Let us now 
comment on (ii) and (iii). Most models studied so far (in particular in [la . Il9|) are based on the 
assumption that the pointer and bath are initially statistically independent. Taking instead the 
whole apparatus to be initially in a local thermal equilibrium seems more realistic. The Markov 
approximation mentioned in (iii) consists in neglecting memory effects for the (reduced) object- 



pointer dynamics. It assumes the decoherence time tA e c to be larger than the bath correlation time, 



a condition not satisfied in some experiments 



12, 



13, 



14 ]. Since decoherence for macroscopic and 



2l|, this approximation is clearly 



even mesoscopic superpositions is faster than bath relaxation 
unjustified if such superpositions arise during the object-pointer evolution. 

We shall assume a certain ordering of time scales. One of these, denoted by T$, is the charac- 
teristic time for the evolution of the measured observable S under the Hamiltonian of the object. 
A second (classical) time scale Tp characterizes significant changes in position of the pointer under 
its proper Hamiltonian (i.e., in the absence of coupling with the object). The initial temperature 
T of the apparatus sets a time scale h/3 = H/(kBT), referred to below as the thermal time {ks is 
the Boltzmann constant). Finally, the object and pointer are put in contact during a time t- m t. For 
a macroscopic pointer, the limit h(3 <C T-p seems difficult to avoid. Similarly, the decoherence and 
object-pointer interaction times t& ec an d iint are small compared with T-p. During an ideal mea- 
surement, the measured observable S may change but weakly under the full (object + apparatus) 
Hamiltonian H, i.e., e < SeT lt ' ~ S for < t < imt^dcc- Only under this condition can an 
eigenstate of the measured observable S be left almost unchanged by the measurement. One has to 
require that (i) the object-pointer interaction Hamiltonian producing the entanglement commutes 
with S (see [3j, |4J, [5j) and (ii) T$ be much larger than t; nt and idee- It is thus legitimate to assume 

tint, tdec *C Tg , tint, *dec, hf3 <^T-p . (2) 

As far as we are aware, this separation of time scales in ideal measurements has not been fully 
exploited in previous works except in Ref. [la]. Unlike in the latter reference, given (J2j) we shall 
not need a further hypothesis on the bath correlation time Tg and its relation with t; n t and tde C . 



A further key input in what follows is the quantum central limit theorem (QCLT) 22J, [23 ] 
which implies Gaussian statistics (Wick theorem) for the bath coupling agent in the pointer-bath 
interaction. This will allow us to study a broad class of pointers and baths, following the approach 
of Ref. 



211 ] . The harmonic oscillator bath linearly coupled to the pointer [24|, l25(] is one member 
of this class, but more general (non-harmonic) baths as well as nonlinear couplings in the position 
X of the pointer will be also considered. It turns out that the decoherence time t^ec may be 
considerably reduced by allowing such non-linear couplings. 

The paper is organized as follows. The model and its different time scales are introduced in 
Sec. HH We discuss the separation of the time scales ([2]) and the ensuing simplification of the object- 
pointer dynamics in Sec. IIIII Sec. [IV] contains a separate study of the two dynamical processes 
producing the entanglement of the object with the pointer and the loss of coherences between 



well-separated pointer readings. That section is pedagogical in character and serves to fix the 
notation and to introduce the relevant time scales; readers familiar with the theory of quantum 
measurement might want to skip the section save for Sees. IIVCI and IIVD1 Our principal results 
are presented in Sec. [V] discussed in Sec. IVIl and finally derived in Sec. IVIIl Our conclusions are 
drawn in Sec. IVIIH Appendices [X] and [B] are devoted to an example for a measurement apparatus 
and to a technical derivation of an approximation for the pointer-bath thermal state. We review in 
Appendix [C] the general properties of two-point correlation functions used in Sees. [VTIVIIl Finally, 
we discuss the QCLT and its consequences (Wick theorem for the bath correlation functions) in 



Appendix [Dj Let us point out that a short report of our results can be found in 26 1. 

Before going on, some remarks may be permitted to put our paper in perspective. Remaining 
within the frame of quantum mechanics and its probabilistic interpretation, we are concerned with 
unitary evolution of the composite system object + apparatus. We discard information about the 
(dynamics of) the microscopic degrees of freedom of the apparatus ( "bath" ) and their entanglement 
with pointer and object by tracing out the bath (see Sees. IIVBI and [3|, |g, |27|]). We so obtain a 
reduced object-pointer density operator with an irreversible evolution. Pointer and object end 
up in the mixed state ([1]) wherein the different pointer states \tp4>) correspond to macroscopically 
distinguishable positions. Such states have quantum uncertainties in position and momentum 
much smaller than the scales of macroscopic readings. Therefore, the irreversible process ([T]) leaves 
a pointer position revealing an eigenvalue of the measured object variable. Over many runs of 
the measurement, the outcome \ip-p) arises with probability |c s | 2 . Similar behavior arises for all 
processes where initially microscopic fluctuations evolve towards macroscopically distinct outcomes. 



A nice example is provided by super fluorescence where 



light 



pulses with substantial shot-to-shot 



fluctuations grow from initial quantum uncertainties [28|, |29|, |30(] . — Let us also recall that quantum 
mechanics is not compatible with the idea that the specific outcome of a single run is predetermined 
by some unknown but "real" property of the object (such a property being independent of the 
measurement apparatus) [311 ] . Competing hidden- variable theories which indulge in such more 



32 



33] 



"intuitive" notions of reality than quantum mechanics have been experimentally falsified 

in ever larger classes, most recently even in non-local variants |34j], while the respective quantum 

predictions were invariably confirmed. 



II. MODEL 



As many authors [3J, \4, \M, la, LL2I> lla uM we consider a three-partite system: the object of 
measurement is some microscopic system (5); a single-degree-of- freedom macroscopic pointer (V) 
will allow readouts; finally, a bath {B) with many (N S> 1) degrees of freedom serves to decohere 
distinct pointer readings. We shall have to deal with the following dynamical variables: for S, the 
observable S to be measured; for V the position X and momentum P; and for B, a certain coupling 
agent B given by a sum of N self-adjoint operators B v acting on single degrees of freedom of the 
bath. The pointer is coupled to S and B via the Hamiltonians 

N 

H vs = eSP , H P B = BX a , B = A" 1 / 2 ^ B v (3) 

v=\ 

where e is a coupling constant and a a positive integer. The object-pointer coupling Hps is chosen 

so as to (i) not change the measured observable S (i.e., [H-ps,S] = 0); (ii) be capable of shifting 

the pointer position by an amount proportional to S, in such a way that each eigenvalue s of S 

becomes tied up with a specific pointer reading; (iii) be a strong coupling (e is large), so that 

different eigenvalues s ^ s' eventually become associated with pointer readings separated by large 

distances. The pointer-bath interaction Hpg is chosen for the most efficient decoherence of distinct 

pointer positions |2ll |. Depending on the value of a, nonlinear (a > 1) as well as linear (a = 1) 



couplings will be considered 35f] . The additivity of the bath coupling agent B in single-degree-of- 
freedom contributions B u having zero mean and positive variance with respect to the bath thermal 
state will allow us to invoke the quantum central limit theorem when taking the limit N — > oo. 
The factor N~ 1 ' 2 in front of the sum in ([3]) is introduced for convergence purposes (the same 
scaling with N is familiar to the classical CLT); note that the pointer-bath coupling constants are 
incorporated within the operators B v . 

The free evolutions of S, V, and B are generated by the respective Hamiltonians Hs, H-p, and 
Hq. We do not have to specify Hs- The pointer Hamiltonian Hp = P 2 /2M + V(X) has a potential 
V(x) with a local minimum at x = 0, so that V'(0) = and V"(0) > 0. The bath Hamiltonian 
His is like B a sum of Hamiltonians acting on single degrees of freedom, Hq = Y^, u Hb,v We thus 
disregard couplings between different degrees of freedom of the bath. The Hamiltonian of the full 
system S + V + B is H = Hs + Hp + H& + Hps + H-VB- A n example of a physical system realizing 
the apparatus V + B is given in Appendix [A] 

We now proceed to describing the initial states allowed for. It is appropriate to require initial 
statistical independence between object and apparatus. The initial density operator ps of the 



object may represent a pure or a mixed state. Two types of initial conditions for the apparatus 
will be considered. The first one, to be referred to as partial equilibrium, is a product state in which 
V has some density operator pp and B is at thermal equilibrium with the Gibbs density operator 
Pb = ^B ex P( — 0Hb)> wherein /? = (fcgT) -1 is the inverse temperature. For this first initial state 
all three subsystems are statistically independent. In the second (more realistic) initial state, the 
apparatus is in thermal equilibrium according to the density operator p^g = Z~g e -Pi H v+ H B+ H vB) _ 
The two initial states of S + V + B are 

p(0) = ps <S> p-p ® Pb partial-equilibrium apparatus (4a) 

p(0) = ps (g) ppg equilibrium apparatus . (4b) 

We further specify the partial-equilibrium state (|4ap by requiring that the probability density 
(x\pp\x) to find the pointer at position x has a single peak of width Ax = A centered at x = 0. 
A momentum uncertainty Ap = 2ttH/\ defines a second length scale A. A macroscopic pointer has 
both A and A negligibly small against any macroscopic read-out scale A c i ass , 

A < 4vrA < A class , (5) 

where the first inequality is the uncertainty principle. We shall also require that 

s-|" m '^ (6) 

which means that the state p-p is not highly squeezed in momentum or in position. As a concrete 
example we may consider a Gaussian pointer density matrix 

(x\pv\x') = -=L= e -(*+*') 2 /(8A 2 ) e -2^( X - X 'f/^ . ( 7 ) 

\/27rA 2 

If V is initially in a pure state then Xxp pp, = J dx dx'(x\p-p\x') 2 = 1, which implies that this state 

has the minimum uncertainty product AxAp = H/2, i.e., A = 4-7rA. 

The Gaussian density ((7|) also arises if V is in a Gibbs state p-p = Z^, e~^ Hv provided that 

the potential V(x) is confining and (3 is small enough. To see this, we note that the pointer 

observables X and P evolve noticeably under the Hamiltonian H-p on a classical time scale Tp, 

which is much larger than all other (quantum) time scales in the model. In particular, Tp is much 

larger than the thermal time, Tp 3> Hf3. As a result, the matrix elements (x\py? \x'} of p-p can be 

approximated by Zp V^^^'))/^- 27 ^-^) 2 /^, wherein A th = 2-Kh((3/M) 1 / 2 is the thermal 

de Broglie wavelength. The reader may recognize in this expression the short-time behavior of 

the quantum propagator (x\e~ ltHv / \x') for t = —ih(3 (see e.g. [3G]). Since the potential V(x) 



has a local minimum at x = 0, it can be approximated near the origin by a quadratic potential, 
V{x) ~ V(0) + x 2 V"(0)/2. Therefore, for small x and x', (x\pj, r\x') has the Gaussian form ([7]) 
with A = A t h = (/3V'(0)) -1 ' 2 and A ~ Ath- It is important to bear in mind the separation of 
length scales Ath <C Ath *C A c i ass . Inasmuch as the pointer classical time scale Tp may be defined 
as T-p = (M /V" (0)) 1 ' 2 , the fact that Ath is much smaller than A t h is equivalent to Tp 3> H@. To 
fix ideas, for Tp = 1 s, M = 1 g, A c i ass = 1 cm, and a temperature of 1 K the above-mentioned 

length scales differ by more than eight orders of magnitude. Hence © and © are well satisfied if 

(eq) 
p v = p), . 

All of these illustrations, including the Gaussian (|7|), are meant to give an intuitive picture. 
What we shall need in actual fact is the quasi-classical nature of the pointer initial state, as 
implied by ([5]) and ([6]), together with the single-peak character of the initial density of pointer 
positions. 

Let us point out an essential difference between our model and the interacting spin model 
of [19l | . Unlike in this reference, 5 is strongly coupled to a single degree of freedom (the pointer 
V) of the apparatus, e.g. with its total momentum P in a given direction (see Appendix |A|) . The 
coupling of S with the other apparatus degrees of freedom (the bath B, for us) is assumed to be 
much weaker and can therefore be neglected, as will be seen in Sec. IVICI Given the separation of 
time scales ^ and our choice of a quasi-classical pointer initial state, the pointer Hamiltonian Hp 
only plays a role in providing an amplification mechanism, as we shall see in Sees. IIIIBI and llVDI 
Hence allowing V to have two or three degrees of freedom, instead of one, would make the notation 
more cumbersome without changing significantly the results. Our results below should also remain 
valid if the bath consists of interacting degrees of freedom (like in a spin chain) provided that the 
(spin-spin) correlations (B^B U ) in the bath thermal state decay more rapidly than l/\p — v\ as 



231 ] . Decoherence 



|/i — v\ — ► oo. In fact, the validity of the QCLT can be extended in this context 

via coupling with a bath of interacting spins and random matrix models for the coupling and bath 

have been considered in 



19, 



37 



We shall study the dynamics of the reduced state of S + V (object and pointer). That state is 
defined by a density operator pps{t) obtained by tracing out the bath degrees of freedom in the 
state of S + V + B, 

pp S (t) = tr t3 (e- itH / h p(0)J tH / h ) . (8) 



Here and in what follows, tr^ refers to the partial trace over the Hilbert space of 
When tracing out the bath we admit the inability to acquire information about it 



27] 



S,V or B. 



III. SEPARATION OF TIME SCALES 

A. Time scales of object, pointer, and bath 

Let us denote by S(t) the time-evolved observable S in the absence of the coupling Hps, i-e., 
for the dynamics implemented by the "free Hamiltonian" Hs- Similarly, let X(t) and B{t) be the 
time-evolved observables X and B when both couplings H-ps an d Hpg are turned off: namely, 

dj(t) = e itH i /h Oj e~ itH J /h , Oj = S, X or B , j = S,VoiB. (9) 

One may associate with the time evolution of X(t), B(t), and S(t) four distinct time scales. The 
time scale Tp = (Af/V'^O)) 1 ' 2 has been already introduced in Sec. [Ill it is the time scale for 
significant evolution of X(t) (or, equivalently, of P(t) = MdX/dt) when the pointer is in the 
initial state pp. The Gaussian form (|7|) for pp and the no-squeezing condition ([6]) make sure that 
Tp is indeed a classical time. 

The bath correlation time Tq is defined with the help of the n-point correlation functions 

h n (t 1 ,...,t n )=tr B (B(t 1 )...B(t n )p% q) ) . (10) 

For simplicity we assume 

tr B (Bp {c B q) ) = . (11) 

Since the bath has infinitely many degrees of freedom, h n (ti, . . . , t n ) decays to zero as \t m — ti\ goes 
to infinity. We define Tq (respectively t B ) as the largest (smallest) time constant characterizing 



22j that for a bath coupling agent B and 



the variations of h n . It follows from the QCLT of Ref. 

Hamiltonian H B which are sums of N independent contributions coming from single degrees of 

freedom, the n-point functions (|1U|) satisfy the bosonic Wick theorem in the limit N 3> 1. This 

means that h n vanishes if n is odd and is given if n is even by sums of products of two-point 

functions, 

h n (t 1 ,...,t n )= Y^ h 2 (t il ,t jl )...h 2 (t in/2 ,t jn/2 ) . (12) 

pairing of {1, . . . , n} 

That manifestation of the QCLT amounts to Gaussian statistics for the bath correlation functions. 
It follows from (|12p that Tg (tg) can be defined more simply as the largest (smallest) time scale 
associated with the variations of ^2(^1^2) = ^2(^1 — £2) as function of t = t\ — t 2 - More precisely, 
h 2 (t) — whenever \t\ 3> Tg and h 2 (t) ~ h 2 (0) whenever \t\ <C t B . Note that with B in thermal 
equilibrium, the thermal time h(3 figures among the decay rates of h 2 and thus ts < h(5 < Tg. 



The time scale T$ is defined in analogy to £g, so as to signal significant variation of the object 
n-point functions trs(S(ti) . . . S(t n )ps). Let us stress that T$ can be larger than the typical inverse 
Bohr frequency h/\E — E'\ of S (here E and E' are two eigenvalues of H$). For instance, if S is 
(or commutes with) the energy H$, then T$ = 00. 

B. Simplified dynamics and initial state 

We assume that the object and pointer observables S(t), X(t), and P(t) do not evolve noticeably 
under the "free" Hamiltonian Hg+Hp during the time span of the measurement, so that tint) ^dec ^ 
T$,Tp. We show now that thanks to this separation of time scales, the impact of H$ and Hp on 
the dynamics can be fully accounted for at times t <C Tg,T-p by modifying the initial states @j 
according to 

p(0) — ► e -it(-ffs+-ffp)A p (0) e it (Hs+H v )/h _ ( 13 ) 

With that slippage of the initial condition accounted for, one makes a small error by otherwise 
dropping H$ and H-p from the total Hamiltonian H in the object-pointer state ([8]). 

Actually, for times t short compared with T$ and Tp, the full evolution operator can be ap- 
proximated by 

e -itH/h ^ e -it(H B +H vs +H VB )/h e -it(H s +H v )/h ^ \t\<^T s ,Tp. (14) 

To justify that simplification we express this evolution operator in the interaction picture with 
respect to Hq = H$ + Hp as 

jtH /n e -itH/h = T expj-^ I dr (H B + eS{r)P{r) + BX{r) a ) \ ; (15) 

here T denotes the time ordering and S(t), P(t), and X(t) are given by Q. Note that for 
|i| <C Tg,Tp these operators are almost constant in time between r = and t = t and may thus be 
replaced in (fT5"|) by S, P, and X. In other words, the right-hand side of (fT5j) can be approximated 
by exp{— H(Hb + Hps + Hpg) / h} , whereupon (fT4"l) is obtained by taking the adjoint and by setting 
£-» -t. 

More specific remarks are in order for each of our two initial states @. We first comment on the 
partial-equilibrium (|4ap . Due to its assumed quasi-classical nature, the pointer state p-p is weakly 
modified by the substitution (I13h in the range of time under study. Deferring the justification of 
that statement to Appendix [B] we shall use e -1 *^ 73 ' h pp e ltHv ' h ~ pp for t ^C T-p. As regards the 



10 




(a) t = 



(b) t S3 tent 



FIG. 1: Wigner function Wp(x,p;t) of the pointer reduced state pp(i); 5 is a spin one-half with two 
eigenvalues ±<5s/2 and (±|ps|±) = 1/2. Note the absence of ripples: the reduced state is a mixture, not a 
superposition of single-peak states. In the horizontal axis, position and momentum are measured in units 
of A/2 and Ap; in the vertical axis, units are such that W-p(x,p;t) has maximum value 1. 



object of measurement S, the free time evolution of the object in 



p° s (t) = e- itH ^ h Ps e itH ^ h 



(16) 



cannot be neglected, even for t <C T$. (For instance, for S = H$ one has T$ = oo and (s\pg(t)\s') = 
e -it(s-s )/ h ( s \p s \ s '} ig n ot close to (s\ps\s') for all finite times t if s ^ s'). However, it suffices for 
our purposes to notice that the diagonal elements (s\p^(t)\s) ~ (s|ps|s) remain nearly unaffected 
by the free evolution when t <C T$. Actually, ti(S(t)ps) = J2 S s ( s \Ps(t)\ s ) nas to approximate 
tr(Sps) = ^2 S s(s\ps\s) in this limit by definition of T$. 

When allowing the apparatus to start out from thermal equilibrium according to (|4b[) . we 
shall take advantage of the "high-temperature" condition h/3 <C T-p discussed in the Introduction. 
We argue in Appendix [B] that under this condition (which to violate for a macroscopic pointer 
would be a nearly impossible task) and for a weak enough pointer-bath coupling satisfying r] t h = 
^-2(0) 1 ' 2 A" h /3 < 1, the Gibbs state of the apparatus can be approximated by 



„(eq) 



ZvB 



e -(3H r /2 e -(3(H B +H VB ) e -(3H v /2 



H(3<.T V . 



(17) 



Moreover, e" 1 "^/ V£b e lMp / R ~ p™ as long as t, hfi <C Tp. 



PVB e 
In conclusion, for both initial states (Ha)) and (14b j) . the substitution (fj~3l) amounts to replacing 

PS by p%(t) in the object-pointer initial state. 



11 



IV. ENTANGLEMENT AND DECOHERENCE SEPARATED 
A. Entanglement of object and pointer 

Before studying the dynamics generated by the total Hamiltonian H, it is instructive to discuss 
what happens if we discard the free dynamics of S, V, and B as well as the pointer-bath interac- 
tion. For the initial state (jjaj), the bath can then be ignored and the object-pointer entanglement 
produced by the interaction H-ps = £ SP becomes particularly easy to describe. Recalling that P 
is the generator of space translations we have e leSPt ' h \s,x) = \s,x — tes), where \s,x) is the joint 
eigenstate of S and X with eigenvalues s and x, normalized as (s,x\s',x') = 5 ss '5(x — x'). Hence 
an initial product state of S + V becomes at time t 

pvs {t) = e- UH ™/ h p s ®p T j tH rs/K 

= Y,( s \PsW)\s)(s'\ J dxdx' (x s (t)\p v \x' s ,(t))\x)(x'\ (18) 

s,s' 
with 

x s (t) = x — tes , x' s ,(t) = x — tes' , (19) 

and, for the Gaussian initial state (|7|), 

{ Xs {t)\p V \x'At)) = -7=L= e -(^'-t,( S+ s'))V(8A 2 ) e -2^ (a: -x'- te ( S - S ')) 2 A 2 . (20) 

v27rA 2 
It is now well to put forth a specification: throughout the present paper we assume for simplicity 
that S has a discrete and non- degenerate spectrum. Moreover, if the Hilbert space of S has infinite 
dimension we restrict ourselves to initial states of the object satisfying {s\ps\s'} = if s and s' 
belong to a part of the spectrum containing arbitrarily close eigenvalues (near an accumulation 
point). 

In the state (|18|) . the diagonal (s = s') matrix elements of the object state ps are multiplied 
by the pointer density matrix p-p shifted by tes in position space, as given by (|20p for s = s'. The 
interaction has thus tied up each eigenstate \s) of S with a pointer state which has position x ~ tes 
with uncertainty A and momentum p ~ with uncertainty 2ttH/\. In position representation, each 
of these pointer states has a peak at x = tes. The different peaks are separated by at least by the 
distance te5s, where Ss is the minimum of \s — s'\ over all pairs (s, s') of non-degenerate eigenvalues 
such that {s\ps\s') ^ 0. In order to be able to infer the value of s from the position of the pointer, 
one must wait until all peaks are well resolved. That resolvability begins at the entanglement time 

tent = -T • (21) 

eos 
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At that time, the reduced pointer density operator pp(t) = trs(f>vs(t)) nas a Wigner function as 
represented in Fig. [2(b). Much later yet, the separation between the peaks reaches a macroscopic 
value A c i ass at the time 

'•class = J --?> *ent > V^2) 

eos 
allowing for a "reading" of the result by a classical observer. 

The entanglement in the state (118p comes from the off-diagonal (s ^ s') contributions in p$. 
Due to the peak structure of the pointer matrix elements (f20|) . for fixed s 7^ s', \(s, x\pps(t)\s', x')\ 
reaches its maximal value when x = est and x' = est. For those values of x and x', 

{s,x = ets\ Pvs {t)\s',x' = ets') = {s\ Ps \s') {0\p v \0) (23) 

is time-independent and proportional to (s\ps\s'). Hence all coherences between different eigen- 
states of S present in the initial state of the object are still alive, no matter how large the time t 
is. At times t > t c i ass , pvsit) resembles a Schrodinger cat state, i.e., has nonzero matrix elements 
between macroscopically distinguishable pointer position eigenstates. For such an object-pointer 
state, no classical probabilistic interpretation is possible: one cannot assign a probability to the 
pointer being located e.g. in the vicinity of x = ets, henceforth implying that S has the value s. In 
a quantum measurement, the entanglement process must be completed by a decoherence process 
suppressing the coherences (j23j) for s ^ s'. 

B. Decoherence and "disentanglement" of object and pointer 

We now turn to the decoherence brought about by the pointer-bath interaction Hpg = BX a , 
momentarily disregarding all other terms in the full Hamiltonian H. As shown in [21] for a similar 
model, a quantum superposition of coherent states of V with well-separated peaks in position 
evolves under H-pB to a statistical mixture of these coherent states. In the situation under study 
here, S and V are entangled, and then decoherence also modifies S. The present subsection 
highlights the fundamental role of this decoherence in a measurement (for more details, see [3j, |5|, |g] ) . 

Assume object and pointer at time to entangled, with p-ps = P-ps gi ven by fT8|) ; the time to 
should be chosen larger than t en t, possibly as large as the classical time scale introduced above, 
tent "C to ~ ^ciass- At time to, the state of S+V+B is p(to) = p^g^PB and the pointer-bath coupling 
Hpb is switched on. To simplify the discussion, let us take for pg a pure state p& = l^oK^oli 
where |^o) = ®i/|Vv) i s a product of iV single-degree-of- freedom wavefunctions. (All arguments 
below can be easily extended to a bath in an initial mixed state like ps = Pg ■) Moreover, let 
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us suppose that (ip u \B u \iJj u ) = and that higher moments {B$) = {ip u \Bl\tp u ) (q = 2,3,...) are 
bounded uniformly in v. The eigenstates \s,x) are entangled at time t > to with the bath states 
\V x (t)) = e -Kt-to)x a B/h^ o y The density operator of S + V + B reads 

p(t) = e ~iH VB (t~t )/h p ent ^ pBe iH VB (t-t )/h 

= Y, f dxdx'{s,x\p^\s',x') \s)(s'\®\x)(x'\®\^ x {t))(^ xl {t)\ . (24) 

s,s' 

We now argue that for x ^ x' , the scalar product {^ x (t)\^ x '(t)) is vanishingly small when the time 
span t — to is larger than a certain decoherence time td ec (x, x'). Due to the additivity ([3]) of B, 

N 

= n (i - (t -'° )V ^; m ' |jm> + o {N -^)) . ( 2 5) 

V = l 

Taking the limit N — » oo for fixed values of t, to, x, and x' we obtain 

<\Mi)|*M*)> = e" D '^') = exp|- ^~^°J 2 )2 +0(JV- 1 /2)| (26) 

t dec (x,x') = (27) 

V ' \x fQ - x^^B 2 ) V ' 

with (B 2 ) = {^o\B 2 \^>o) = -W -1 X^W^I-^lVv}- We have so far retraced the proof of the (classical) 
central limit theorem. 

Taking the partial trace of (|24p over the bath Hilbert space yields 

PTsit) = ^ j dxdx' (s\p s \s')(x s (t )\pv\x' s ,(to))e- Dt ^ \s)(s'\ ® \x)(x'\ . (28) 

s,s' 

Due to the coupling with the bath, each matrix element of ptpg is now multiplied by the scalar 
product (|26p . Let us consider a particular term s ^ s' in the sum in the right-hand side of (|28p . 
To simplify the forthcoming discussion, we assume that a = 1. As follows from the peak structure 
of the pointer coherences (|20p . only the terms satisfying x ~ x s o = toes and x' ~ x s /q = toes' with 
uncertainty A contribute significantly to the integral over x and x'. For those terms, Dt(x,x') = 
D t (x s o,x s io)(l + O(t cnt /to)), see (f2~Tj) . ([26]) . and (p7|) . Therefore, if t — to is large compared with 
tdcc(x s o,x s 'o) and to 3> t en t, the product (x s (to)\p-p\x' s ,(to)) e~ Dt ( x ' x ' is vanishingly small for all 
values of (x,x'). The off-diagonal terms corresponding to s / s' then become negligible in the 
object-pointer state (|28|) . It is worth emphasizing that tdec( x sO, Xs'o) can De niuch smaller than the 
dissipation time scale on which the pointer-bath coupling irreversibly changes the pointer position. 
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We would like to point out that the aforementioned damping of the coherences is related to a 
lack of information about the bath in a more subtle way than what is suggested by the partial 
trace in (|28p . In fact, some partial knowledge of the bath state would not inhibit this decoherence. 
More precisely, in order to obtain some information on the coherences in the full density operator 
(terms proportional to \s)(s'\ with s ^ s' in (JMD) at times t — to 3> td C c( x sO, x s /q), it is necessary to 
perform a measurement on some bath observable Ob satisfying (^f Xs0 (t)\OB\^x a , Q (t)) ^ at such 
time. It can be shown by repeating the arguments yielding to (|26p that such an observable must be 
non-local, i.e., it must act non-trivially on all bath degrees of freedom except for a finite number of 
them. Considering that measuring such an observable is "unrealistic" , everything happens as if the 
s ^ s' terms have disappeared in ([2"3]) . The crucial point is that the object-pointer state is entangled 
by H-pQ with a very large number N of bath variables, so that information about the coherences is 
spread out between these many variables after some time. Macroscopically distinguishable object- 
pointer states are then entangled with bath states which are almost orthogonal in many subspaces of 
the bath Hilbert space. This makes the situation quite different from the entanglement discussed 
in Sec. IIV Al there, the object state was entangled with a single pointer variable x and it was 
implicitly assumed that any pointer observable (in particular, its position X) could be "observed" 



30] (Sec. 22.11) 



, y, m 



at some ultimate stage of the measurement. We refer the reader to 

(Chapter 2), and [6] for related discussions on this very important conceptual point. 

Let us define the decoherence time idee as the largest of the times idec^sOi^s'o) f° r all distinct 
eigenvalues s and s' such that (s\ps\s') / 0. For t — to S> tdeo the object-pointer state has shed all 
terms s ^ s' in the double sum in the density operator (|18j) . 

Pvs{t) ^ ^(s|As|s> |s>(a| ® Pv(t) (29) 

s 

wherein it has been assumed that idee *C t — to <C t&,Ts,Tp (so that Hg, Hs, and Hp can be 
neglected) and t en t *C to and we have set 



Pv 



it) = f dxdx' (x s (to)\p v \x' s (to)}e- {t - to)2/td ^ x ' x ' )2 \x)(x'\ . (30) 



While pps{t) is n °t (and actually never can become) strictly diagonal in the position basis of the 
pointer, the matrix elements of the pointer state (|30p almost vanish if \x — x'\ is larger than either 



the uncertainty A (see (|2"0]) ) or the decoherence length y/2h/(a A a_1 (t — to)y / (-B 2 )) (see the second 
factor inside the integral in (f30fl ). 

It is worth noting that the object-pointer states appearing in the sum over s in (|29l) are product 
states; pps(t) is a statistical mixture of these states with probabilities p s = (s\ps\s). Hence the 
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decoherence disentangles S and V. This implies that in the time regime indicated after (|29p . S 
and V can be given independent states ps(t) and pv(t), 



Ps 



(t) = Y,Ps\s)(s\ , Pv (t) = Y,PsPv(t)- (31) 



The object S is in one of the eigenstates \s) with probability p s , in agreement with von Neumann's 
postulate. The pointer V is in the quasi-classical state Pp(t), with the same probability. 

C. Summary 

Let us sum up the discussion of the two previous subsections about the object-pointer entan- 
glement produced by the interaction Hp$ and the decoherence arising from the coupling with the 
bath H-pb- The dynamics implemented by Hp$ uniquely ties up after the entanglement time t en t 
each eigenvalue s of S with a characteristic pointer position x s (t). Such neighboring pointer posi- 
tions differ by more than the uncertainty A then. Note that arbitrarily close eigenvalues cannot be 
resolved within a finite time: in fact, i e nt tends to be large for an object initially in a superposition 
of eigenstates \s) with closely lying eigenvalues s (e.g. near an accumulation point of the spectrum), 
i.e., for small values of 5s in IJ21J) : this limits in practice the precision of the measurement of S. 
In the absence of any other interaction, after a time i c iass S> tent the initial product state of the 
object and pointer has evolved into a Schrodinger cat state. Nothing irreversible is brought about 
by the dynamics: the entanglement can be as easily undone as done, by applying the Hamiltonian 
H-ps with the parameter reset e — ► — e. 

The dynamics generated by Hpg brings about decoherence. After the decoherence time t^ ec , 
any pair of object-pointer states corresponding to macroscopically distinguishable pointer positions 
are entangled with almost orthogonal bath states. After averaging the object-apparatus state over 
the bath variables, one obtains an object-pointer state pps(t) in which all information about the 
coherences between such states is missing, i.e., all coherences for pairs (s, s') of distinct eigenvalues 
are suppressed. The irretrievable loss of information about the bath goes hand in hand with the 
irreversibility of the object-pointer dynamics. 

An object-pointer state pps{t) describes an accomplished measurement under two conditions: 

(i) All coherences (s,x\pps(t)\s' ,x') corresponding tests' have disappeared, so that S + V is 
in a statistical mixture of separable states like in (|29|) : this occurs at time t 3> idee- 



(ii) The separation between the peaks of the distinguished pointer states p-p(t) reaches a macro- 
scopic value A c i ass ; this occurs at time to ^ ^ciass j see (f22j) . 
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FIG. 2: Sketch of a candidate for the pointer potential. The height Vo of the potential barriers around x = 
and the width W of the potential wall are much larger than ksT and the thermal fluctuation Ath- 

Only for to ^ ^ciass can a classical observer infer a measured value s by looking at the position of 
the pointer. Such a "reading" of the pointer, while still a physical process in principle perturbing 
V, surely cannot blur the distinction of the peaks. Rather, the pointer will behave classically under 
a reading, i.e., it will not noticeably react. 

D. Unstable pointer potentials and amplification 



It is clear from (|22p that condition (ii) can hardly arise unless the object-pointer coupling 
constant e is very large. This is related to the well-known amplification problem in quantum 
measurements 39(. In order to get rid of this unrealistic condition on e, one may consider a 



different situation than that described in Sec. HVAl Let us take a non-confining pointer potential 
V(x) with two potential barriers separated by a distance W, see Fig. [2J The height of these barriers 
is large compared with the thermal energy 1/(3. We now replace the initial states pp = p-p in 
(|4ap and p-pg in (|4bj) by local equilibria within the potential well. (This local equilibrium for the 
apparatus can be achieved by first preparing V in some state localized near x = at time t = —t\, 
with £; larger than the relaxation time but small compared with the tunneling escape time, and 
then letting V interact with B until t = 0.) Our previous statements about the distinct peaks in 
the pointer density produced by the object-pointer interaction remain valid for such initial states. 
The interaction Hps is switched off at some time tint- If tint is larger than W/(eSs), the separation 
between the peaks in the pointer density at time t; n t will be subsequently amplified by the pointer 
dynamics. Assuming also that A <C W <C e5sT-p, A c i ass , one has t- mi w W/(e5s) <C 7p,t c i ass . In 
this situation, the small quantum system S must be able to perturb the pointer strongly enough in 
order to produce in it a "mesoscopic change" (i.e., a distance W between the peaks in its density), 
instead of a macroscopic change as required in Sec. IIVA1 In particular, if p-p is a Gibbs state 
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with position uncertainty A = A t h = (,3V"(0)) _1 ' 2 , this arises when the height Vq ~ W 2 V"(0) 
of the two potential barriers satisfies (3~ l < V < M(e5s) 2 and y < ^dass V "(°) (recall that 
Tp = M/V"(0)). Then condition (ii) of the preceding subsection will be fulfilled after the object- 
pointer interaction has been turned off, at time to ~ T-p. If moreover Vq <C V"(0)(e5sTs) 2 , the 
simplification of the dynamics discussed in Sec. lIIII can be used since the object-pointer interaction 
time satisfies £; n t <C T$,Tp. 

V. SIMULTANEOUS ENTANGLEMENT AND DECOHERENCE 

We now present and discuss the main results of this work, before deriving them in Sec. IVII1 We 
are interested in the object-pointer dynamics when, unlike in the situation just described, S and 
V evolve under the simultaneous action of H-ps arm H-pg. Furthermore, in contrast to Sec. IIVBI 
we do not neglect the bath Hamiltonian H&. 

A. Partial-equilibrium initial state 

Let us first consider the evolution of the initial state (jjaj). Due to both the initial statistical 
independence and our special choice of the interactions, the density matrix of S + V retains at 
"short" times t <^T$,Tp a remarkably simple product structure (see the discussion in Sec. IIII|) . 



{s,x\pp S (t)\s',x') = {s\p° s (t)\s') {x s (t)\pp\x' s ,(t)) exp{-D t (x s (t),x' s ,(t);s,s')-icp t } . (32) 

Here p$(t) is given by (fT6j) . x s (t) = x — tes, x' s ,(t) = x' — tes', and 4>t is a certain real phase 
(depending on t, x, x' , s, and s') which we do not specify here since it is irrelevant for decoherence. 
We shall derive in Sec. IVIIl the decoherence exponent D^, 

D t (x,x';s,s') = ^ f d Tl f d^dx' + nes'r-ix + nesr) x 

(Or' + T 2 es') a -{x + T 2 es) a )h( Tl - r 2 ) (33) 



where h{j\ — t 2 ) = h 2 (Ti, t 2 ) is the bath two-point function defined in (|10p . The first factor in (|32|) 
accounts for free evolution of the object initial state ps, as generated by H$, see (fl"6j) . It is equal 
to p s = (s\ps\s) if s = s' and t < Is. The second factor is nothing but the matrix element (|20l) 
of the shifted pointer initial state. Here, the Hamiltonian Hp does not show up because of our 
assumption t <C Tp and our choice of a quasi-classical initial state pp. Most important is now the 
third factor in (j32|) : it accounts for decoherence, i.e., for the suppression of coherences with respect 
to pointer displacements associated with different eigenvalues s / s'. 
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The exponent Dt has the following properties: 

(a) D t (x, x'; s, s') > for all values of x, x' , s, and s'. 

(b) Dt(x,x';s,s') = initially (for t = 0) for all matrix elements and at all later times for the 
diagonal matrix elements (x = x' and s = s'). 

(c) D t (x s (t),x' s ,(t);s,s') = D- t (x,x';s,s'). 

The non-negativity (a) is a consequence of the fact that the correlation function h(t) and its real 
part 3ft/i(i) are of positive type, i.e., they have nonnegative Fourier transforms h(w) and (Jfth)(u). 
Actually, one may rewrite ([33]) as 

1 f°° c\lj . /** 2 

D t (x,x';s,s') = -- —^h)(io)dT((x' + TesT-(x + re S ) a )e-^ T >0 (34) 

where we have used (3ft/i)(u;) = (3?/i)(— cj) and (3>h)(u}) = —(Qh)(—u). Property (c) is easily 
checked by a change of the time integration variable in (|34|). Let us recall from Sec. IIV Al that 
the dynamics generated by H-ps maps the object-pointer coordinate (x,s) to (x s (t),s) after time 
t and, similarly, (x',s') is mapped to (x' ,(t),s'). Hence one may interpret (c) as the invariance of 
Dt under time reversal, i.e., under t — > — t and the exchange of the initial and final coordinates. 

B. Equilibrium apparatus initial state 

Our result for the initial state (j4b|) looks quite similar to that for the initial state (|4a|). Before 
stating it, let us introduce the effective pointer potential 

V ea (x) = V(x)-h- 1 l0 x 2a , (35) 

wherein 70 is given in terms of the imaginary part of the bath correlation function h(t) by 

,0 
7o = / dr^h(r). (36) 



We write Z-p e g = J dxe~^ Vea ^ x ' for the partition function associated with V e R. It follows from the 
general properties of h(t) that < 70 < h(3h(0)/2 (see Appendix [C]). Considering e.g. a linear 
pointer-bath coupling, 2jq/H is the mean force per unit length exerted by the bath on the pointer. 
Note that V e s(x) is a non-confining potential if V(x) = o(x 2a ) at large distances. For instance, 
if V is a harmonic oscillator (V(x) oc x 2 for all x) and a > 1 then V e g(x) looks like in Fig. [2j 
This means that an initial pointer density localized around x = will tunnel away and eventually 
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spread over the whole real line once the pointer-bath coupling is switched on. In such a case the 
apparatus equilibrium state ppg must be replaced by a local thermal equilibrium (see Sec. lIVD]) , 
This local equilibrium exists under certain conditions on the pointer-bath coupling to be discussed 
below. 

As we shall show in Sec. IVIICl the object-pointer density operator is given at times t <C Tg, T-p 
by 

(s,x\p vs (t)\s',x') = (s\p°s(t)\s') R t (x s (t),x' s ,(t);s,s') exp{-D t (x s (t),x' s ,(t);s,s') -ifa) (37) 

with the same decoherence exponent Dt and phase (ftt as above. The only difference between 
(|37p and the formula (|32p for the partial-equilibrium initial state lies in the replacement of the 
initial pointer density (x\pp\x') by the function R t (x,x';s,s'). For a time t short enough so that 
Dt(x, x'; s, s') < 1, this function is given by the Gibbs-type density 

Rt(x, x'; s, s') ~ R (x, x') = Z~)^ e ^0^(*)+^(*'))/ 2 e -2^'-) 2 /A t \ (38) 

For larger times t (with the proviso t <C T$, T-p), Rt is given by the more complicated integral (|85p 
or, in the special case a = 2, by the formula (|89|) below. Let us only mention here that for a = 1, 
([38[) gives the correct answer up to a phase factor for all times t <C T$,T-p. Interestingly, ([37]) 
entails the following result on the reduced pointer initial state 

(x\tr B (p^)\x')=Ro(x,x'). (39) 

Comparing (|38|) with the expression of (x\p!p \x') at high temperatures given in Sec. [Til we see 
that the coupling between V and B can be fully accounted for by the effective potential (|35p . 
Furthermore, for a linear coupling a = 1 the matrix elements (|38|) can be approximated for small 
x and x' by the Gaussian ([Jj) with an almost unchanged uncertainty in momentum, Ap ~ 2irh/\ t h, 
and a renormalized uncertainty in position A c g > Ath given by A^j = /3V^(0) = /3(V" (0) — 2^/o/h) 
(see Sec. HI]) . 

Our results (|37H39f) rely, in addition to t <C Tp,T$, on two additional hypotheses: (a) the 
separation of time scales H/3 <ti Tp or, equivalently, the separation of length scales Ath *C Ath (see 
Sec. IIHB[) ; (b) a weak enough pointer-bath coupling satisfying 

r/th < l/v 7 ^ if a = 1 

with r ?th = (B 2 ) 1 / 2 A t " h /?. (40) 

r/th Cl if a > 1 

Here (B 2 ) = trg(5 2 /9g ) = /i(0) is the thermal variance of the bath coupling agent. Condition (j40p 
is motivated by the following requirement: The effective potential (j35p must have a local minimum 
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at x = and the height of the potential barriers surrounding the origin must be large compared 
with the thermal energy 1//3. Only under that condition can pointer and bath be prepared in a 
local thermal state in which the pointer reduced state has a single peak at the origin like in Fig. 
QJa). If the coupling H-pB induces an instability in the pointer-bath dynamics, we must replace 
the Gibbs state p^J in (|4b|) by that local thermal state, as explained in Sec. IIVD1 Note that we 



mm 



exclude here pointers being at a critical point of a phase transition considered in 

We first consider the case a = 1. If V(x) = V"(0)x 2 /2, the aforementioned requirement is met 
whenever V^j(O) > 0, i.e., Jq/H < V"(0)/2. This stability condition is well known for a harmonic 
oscillator interacting linearly with a bath of harmonic oscillators [40]. For a potential V(x) which 
is non-quadratic at large distances |x| > W, we must stipulate a bit more, e.g. jo/fr < V"(0)/4, in 
order that the height of the two potential barriers be large compared with 1/(3. Bearing in mind 
that 70 < hf3(B 2 )/2, the latter condition is satisfied under our hypothesis (|40p . Most importantly, 
it implies A t h < A g < y^A^, so that the various length scales are ordered as Ath *C A t h as 
A cff < W. 

Now turning to the case a > 1 we insert V(x) — V"(0)x 2 /2 into (j35[) and find a distance 
between the left and right maxima of V e g(x) equal to W e g = 2(hV"(0)/(2ajo)) 1 /( 2a ~ 2 \ these 
maxima equaling V"(0)(hV"(0)/ , yo) 1 '^ a ^ 1 ' up to a factor of the order of unity. As a result, (|40p 
implies the required stability of V e g(x). According to the discussion of Sec. IIV Dl the object-pointer 
coupling can be switched off at time £i nt ~ W e g/(e5s) fa (hV"(0)/ r yo) 1 ^ 2a ~ 2 \^s)^ 1 . This time 
must be chosen small compared with T$ and T-p and large compared with the entanglement time 
t ent = A th /(e<5s), so as to fulfill © and flUJ). 

By comparing (|32|) and (|37|) we may conclude that the coherences of pj>s(t) for s ^ s' decay 
to zero in same way for the two initial states (JH), at least in the early time regime when these 
coherences are not yet very small. Furthermore, in view of (|38|) the whole discussion of Sec. IIVI 
about the emergence of classically discernible peaks remains qualitatively valid. 

VI. DECOHERENCE TIMES 

Before presenting a derivation of our main results (I32h and (|37|) in the next section, we focus our 
attention to the decoherence factor e~ Dt . It has been stressed in Sec. IIV B1 that the object-pointer 
matrix elements 

p P v f(t; S , S ') = (s,x = ets\p PS (t)\ S ',x' = ets') = {s\p% [t)\s') R t (0, 0; s, s>) e^f^M-^ (41) 
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are of particular importance for decoherence in a quantum measurement. Here Rt(0, 0; s, s') is 
equal to (0|/>p|0) for the initial state (|4"a|) . to Z^ cS for the initial state (|4b|) if Df ca (s,s') < 1, 
and to a more complicated function of s and s' for the initial state (|4bp if D^ oa (s,s') > 1. The 
main difference between (|4"Tj) and (f23j) lies in the presence of the damping factor exp{— Df ea (s, s')} 
given by 



2« ft ft 

Dr k (s,s') = D t (0,0;s,s') = w {s' a - s a ) 2 J d n J dr 2 rf r^n - r 2 ) . 



(42) 



A. How does _D^ ca grow with time? 

The decoherence factor (1421) is positive, vanishes for s = s' (see (a) and (b) in Sec. IV All , and 
satisfies the following properties: 

(d) D\ ca (s, s') is an increasing convex function of time if s a ^ s' a . 

(e) Df 03, (s, — s) = if a is even. 

(f) A0r,x'; s, a') = A P ° ak (s, a') (l + 0((M + |x'|)(ei|s - s'D^ 1 )) for |x|, |x'| <C et|a - s'|. 

Property (d) means that, quite generally, the graph of Df 03, looks qualitatively like in the inset 
in Fig. [3l To establish this result, we take x = x' = in (|34p . differentiate both sides with respect 
to t, and do the time integration by parts to get 

Ot it,- Jo 7I w Jo 

Using the fact that the function (1 — u)° _1 is positive and decreasing between and 1, it is easy 
to show that the integral over u in (|43p is positive for almost all w > 0. Bearing in mind that 
($lh)(u)) > 0, this establishes that dDf CS/ jdt > for t > 0. Hence Df ea is an increasing function 
of t. By a similar argument, d 2 Df ca /dt 2 > and thus Df ea is convex. 

According to property (e), if a is even and the spectrum of S is symmetric with respect to 
s = 0, the coherences (j41|) for s' = —s are not damped. This comes from the symmetry x <-> — x of 



9 npoak/ r\ e a I ia q\2 ,2a f°° duJ (^th)(uj) f j . 

— iJf ( S ) S ) = T2"( S — s ) ai / / du(l — u) sm(ojtu) . (43) 

at n J it u) Jq 



the Hamiltonian H-pg in ([3]) , which allows for the existence of decoherence- free subspaces 4l( . Due 
to these long-living coherences, V + S fails to reach (at least within a time span t <C Tp,Ts) the 
statistical mixture required to be able to give a classical result to the measurement. We exclude 
that case from now on. More precisely, we assume that if a is even then s/s' is not close to —1 
for all pairs (s, s') of eigenvalues such that (s\ps\s') ^ 0, i.e., \s' a — s a \/\s' — s\ a is bounded below 
by a constant c™ in > of the order of unity. With this restriction, for (\x\ + |a/|)/(e£|s — s'\) 
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sufficiently small the error term in property (f) is bounded by (\x\ + |x'|)(ei|s — s'|) _1 times a 
constant independent of x, x', s, s', and t [421 ] . 

We are concerned in this section with determining the time scale td ec ( s i s ') characterizing the 
growth of Df ea (s,s r ) and the corresponding decay of the (s 7^ s')-coherences (|4ip . This time, to 
be called the decoherence time, is defined implicitly as Df™ (s, s') = 1, i.e., 

l~wJ = ^rL dri /o dT2TiT2 W) (44) 

where 

tent(s,s') = — (45) 

e\s — s\ 

is the entanglement time (whose physical interpretation has been illustrated in Sec. lIV A\\ , 77 is the 
(fluctuation of the) initial pointer-bath coupling energy in units of A^T, 

17= (By/ 2 A"^(3(tv(Hl, B p v ®p^)) l/2 , (46) 

c a (s, s') = 1 if a = 1, and 

c a (s,s')= l — J-! if a>\. (47) 



For the initial state (|4b|) . one must set A = A t h in (|45p and (|46p and 77 = 7/ t h must be small 
enough, see (j40j) . By inspection of (f44|) . tdec(s ) 'S / ) depends on the object-pointer and pointer-bath 
coupling constants e and 77 through a single parameter er] l ' a . Recalling that (|32M42p are valid with 
the proviso t <C Ts, Tp, the "free" evolutions of S and X must be slow compared to td ec (s, s'), i.e., 

t dec (s,s')«£T s ,T v , s/s' . (48) 

For given s ^ s', if td ec (s, s') > t cnt (s, s') then at time t 3> t^ cc (s, s') the peaks at (x, x') = (ets, ets') 
of the pointer coherences in (|32p and (|37p (second factors on the right-hand sides) are flatten down 
by decoherence (third factors), so that (s, x\p-ps{t)\ s ' , x') ~ for all values of (x, x'). This statement 
follows from a similar argument as in Sec. IIVBI and from property (f ) (see the beginning of this 
section). It is worth emphasizing that if, unlike in the situation just described, £dec( s > s ') is smaller 
than t en t(s,s') then the coherence (s,x\p-ps(t)\s',x') may still be large at time td ec (s,s') for some 
(x, x') ~ (ets, ets') with uncertainty A. In such a case the decoherence time must be defined as the 
time t at which the minimum of Dt(x, x'; s, s') over all values of (x, x') is equal to 1. We postpone 
to a separate work the determination of that decoherence time. 
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The decoherence time idee of the measurement is the largest of the times t^ ec (s, s') for all pairs 
of distinct eigenvalues (s,s') such that (s\ps\s'} ^ (with the proviso idee > *ent = A/(e6s) in 
light of the discussion in the preceding paragraph) . This amounts to replacing | s' a — s a \ in (|42[) by 
its minimum value over all such pairs (s,s') (recall that Df 03, is an increasing function of time). 
For a = 1, this minimum value is by definition equal to 5s (Sec. IIV Ap ; for a > 2, it depends on 
the spectrum of S in a more subtle way |43J|. At times t 3> t^ ec , the object-pointer state p-ps{t) is 
very close to the separable state (|29p . In other words, S and V are in the statistical mixture (|31|) 
with the probabilities p s = (s\ps\s) and with pointer states Pp(t) given by 

(x\pf>(t)\xf) =Rt(x s (t),x' s (t);s,s) exp{-D^ t (x,x';s,s)-i4> t } (49) 



with Rt(x , x 1 ; s , s) equal to {x\p-p\x') for the initial state (|4"a|) and to the Gibbs-like density ([38|) 
for the initial state (J4b]) when D_t( x , x ']s,s) < 1. We have used in (|49p the time-invariance 
property (c), see Sec. IV Al The initial superpositions of eigenstates \s) have disappeared by indirect 
decoherence via the pointer. The pointer is in a statistical mixture of quasi-classical states having 
densities localized around x = tes with uncertainty A. The essence of quantum measurements 
lies in this loss of coherences: for indeed, as already pointed out in Sec. IIV CI it is only when all 
object-pointer coherences for s ^ s' are vanishingly small that a classical probability can be given 
for the result of the measurement. 

The pointer matrix elements (|49p are also damped by decoherence via the last exponential factor 
in (|49p . One can show, however, that for relevant values of x and x' satisfying \x — ets\ < A and 
\x' — ets\ < A, the corresponding damping time is much larger than ideo at least in the two limiting 
regimes idee <C t& and tdec S> Tg studied below. The special case a = 1 will be discussed in the 
next subsection. 

It is worthwhile mentioning here that one should expect that idee *C t c lass> save f° r extremely 
large object-pointer coupling constants e. Object and pointer are then never in a Schrodinger cat 
state as in (fTHI) . because decoherence subdues linear superpositions (via the third factors in (|32|) and 
(|37D ) faster than entanglement between V and S can produce them (second factors in (|32p and (I37p ). 
Due to the simultaneous action of H-ps an d H-pg, the whole measurement process directly produces 
the mixture of macroscopically distinct pointer states p4>(t), without allowing for the intermediate 
appearance of macroscopic superpositions. This is one of the central results of the present paper. In 
the situation described in Sees . IIV Dl and |VB"1 i.e., if the (effective) pointer potential is unstable and 
the object-pointer interaction is switched off at time tmt ~ W/(e6s) <C ?V,t c i ass , even mesoscopic 
superpositions do not appear at any stage of the measurement when idee *C tint- 
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FIG. 3: Decoherence against entanglement times in units of Tg = hj3 for the harmonic oscillator bath 
considered in Sec. IVIDI with r\ = 10 _1 , c = 1, and wd — 5 (log-log scale). Solid curves: exact results 
for (a,m) = (1,5), (1,3), (1,1), (2,5), (2,3), and (2,1) (from left to right). Broken lines: approximate 
expressions for Td cc <C w~^ (dashed curves) and Td cc ^> 1 (dotted curves), see text. Inset: decoherence 
exponent Df ca as a function of r = t/Tg for (a, m) = (1, 3). 

No assumption whatsoever was made on the bath correlation time Tq to establish (|32M42p . Our 
results therefore go beyond the so-called Markovian limit which would require Tq <C idee- This 
is an important point, since for sufficiently large e decoherence will take place within the "non- 
Mar kovian" regime t < Tg. Explicit asymptotical results for idee can now be drawn from the 
foregoing expressions for both idee ^ ^B and t& cc 3> Tq. 



B. Interaction dominated regime i^cc *C ig 



In the (non-Mar kovian) regime idee <C ts, the dynamics is dominated by the interactions H-pg 
and H-pQ. For t <C tg, one may approximate h{r) in (|42p by the thermal variance h(0) = {B 2 ) of 
the bath coupling agent. This yields 

2«+2 



Dr k ( S , s ') 



t 



, '•deev^) ^ ) 



(interaction-dominated regime) 



t ^ M .(^±l>Y'(^il)^ k , 



c a (s,s') 



hPrj 1 /" 



(50) 
(51) 



where we have used the entanglement time and dimensionless parameters (145M7p . That result 
makes sense with the proviso £dec(s, «') <C tg,Ts,Tp. The fact that Df c& oc t 2a+2 could be 
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expected from (|26p : indeed, a contribution of t 2 (x' a — x a ) 2 to the decoherence exponent was found 
for H-ps = and fixed x, x'; recalling that for Hp$ given by ([3D the positions of the peaks grow 
proportionally with time, we recover the above-mentioned power law. It should also be noted that 
idee depends on the bath through the single parameter r]. This comes from the fact that the bath 



dynamics can be ignored when td ec <C tg |44i |. 



Invoking (|42p and (3ft/i)(w) > 0, it is easy to show that (|50p gives an upper bound on Dj? ea (s, s') 
for all times t > 0. By property (d) in Sec. IVI Al in the regime td ec ^ ^B the decoherence time 
must be larger than the right-hand side of the asymptotic formula (|5ip . 

It has been stressed above that the interpretation of idee as the decoherence time of the mea- 
surement relies on the assumption t en t — ^dec • We now argue that this condition is fulfilled if the 
pointer-bath coupling energy is of the order or smaller than ksT (i.e., rj < 1) and s' / s is not very 
close to unity. In fact, under these assumptions one has even t en t(s,s') <C t^ cc (s,s'). This follows 
from (f5"Tj) . the consistency condition tdcc(s,s') <C tg and the inequality tg < hf3, which imply 
td ec (s, s') <C hf3 and thus t ent (s, s') <C hf3. In contrast, if a > 2 and \s'/s — 1| <C 1 one sees by in- 
spection of (|47j) that c a (s, s') ~ a\l — s' /s\ 1 ~ a S> 1. Hence the first factor on the right-hand side of 
(J5T|) is small and one may have £ e nt(s, s') > t<i ec (s, s'). This corresponds to an initial superposition 
of eigenstates |s) with closely lying eigenvalues, as discussed in Sec. IIV Cl Similarly, one may have 
t e nt(s,s') > t^ ec (s,s') for a strong pointer-bath coupling, i.e., for r] S> 1. 

It is worthwhile comparing the strength of decoherence for different values of the exponent a in 
the coupling Hamiltonian H-pg, keeping its magnitude r]//3 constant. We find that idee is smaller 
in the nonlinear case a > 1 in comparison with the linear case a = 1 by a factor t dec /t de ~ of 
the order of c a (t en t/t\ e ~ ) ^ a ^ 1 ''( a+1 > <^ 1. Interestingly, a linear pointer-bath coupling is 

much less efficient than a nonlinear one in suppressing the coherences ()4ip for s ^ s'. This has the 
following important consequence: for a pointer-bath coupling of the form H-pb = B f(X) with f(x) 
a smooth real function, a dipole-like approximation consisting in linearizing f(x) may lead to an 
over-estimation of the decoherence time tdec even if /"(0)A//'(0) is small. Actually, the quadratic 
coupling Bf"(0)X 2 /2 gives a smaller decoherence time than the linear coupling Bf'(0)X when 
tent/h/3 < rj-Hf'^A/f (0)) 2 with rj = /'(0)( J B 2 ) 1 /2 A /?. 

We can now give an explicit condition ensuring that tdec is smaller than the time ti n t ~ W/(e5s) 
needed by object-pointer entanglement to produce superpositions of pointer positions separated by 
the mesoscopic length W, A <C W <C A c i ass : t en t must be large compared with (A/W) a+1 hf3 / (r]c a ) ■ 
In this limit decoherence is so fast that these mesoscopic superpositions do not appear at any 
moment during the measurement. 
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As pointed out in Sec. lVI Al it is appropriate to demonstrate that the decay of the pointer matrix 
elements (|49|) remains negligible for times t until well after the disappearance of the off-diagonal 
(s ^ s') terms in p-pg(t). Due to the peak structure of the pointer density (first factor in the right- 
hand side of (fl9|) ). the relevant values of x,x' are such that |x s (t)|, |x' s (t)| < A. Given t ent <C idee: 
such x and x' are separated by a distance \x — x'\ < 2 A much smaller than the interpeak distance 
t^ cc e5s relevant for the decay of the (s, s' = s + 5s) matrix elements of pvs(t). We restrict ourself 
to the case a = 1 and consider the limit t <C tg < h(5. Setting s = s' , inverting the sign of the time 
t, and replacing h(n - r 2 ) by (B 2 ) in ([33]), one finds that -D_ t (x, x'\ s, s) = {B 2 )t 2 (x - x') 2 /{2h 2 ). 
Note that this decoherence exponent is the same as in (130p . If 77 < 1 and x, x' , and t are in the 
range mentioned above then L>_ 4 (x,x'; s, s) < 2n 2 t 2 (h(3)~ 2 <C 1. Hence the decoherence caused by 
the pointer-bath coupling has a small effect on the pointer states Pj>{t) up to times t < £dec( s )j the 
decoherence factor in (|49p being still close to unity. So indeed, the bath does away with the "off- 
diagonal" (s 7^ s') object-pointer matrix elements before the "diagonal" ones change noticeably. 

C. Markov regime idee ^> Tg 

When idee 3> tg the off-diagonal matrix elements (|4ip have no time to decay between t = and 
tg. Decoherence may then take place within the so-called Markov regime t S> T&, also known in 



the mathematical literature as the singular-coupling limit [4=51 . I4a | . Note that under our condition 
idee "C Ts, Tp it is not appropriate to use a rotating- wave approximation. Decoherence is governed 
in the Markov regime by the small- frequency behaviors of the Fourier transforms (dth)(u) and 
(Qh)(uj) of the real and imaginary parts of the bath correlator h(t). We shall make use of a few 
properties of these Fourier transforms, which are explained in more detail in Appendix O We 
assume that (3 : /i)(o;) ~ — ijto" 1 for u> <C Tg , 7 being a positive constant. Bearing in mind that 
(^sh)(u>) is an odd function of a; and must be regular enough (i.e., admit differentials of sufficiently 
high orders) in such a way that Qh(t) decays rapidly to zero as t — > ±00, we take m to be a positive 
odd integer. By analogy with the case of a bath of harmonic oscillators linearly coupled to V, we 
speak of Ohmic damping when m = 1 and of super- Ohmic damping when m > 1 [24J, [25J] . The 
behavior of (3ft/i)(u;) at small frequencies can be deduced from that of (Qh)(oj) thanks to the Kubo- 
Martin-Schwinger (KMS) relation (|C6p , Such a relation holds because the average in the correlation 



function h(t) is taken with respect to a bath Gibbs state 



47fl. It implies (3tt)(w) ~ 2juj m - 1 /(hf3). 



Let us first discuss the super-Ohmic case m > 3. The frequency integral in (|43p can be rewritten 
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after an integration by parts as 

f°° doo (Wi)(u) f 1 . /n ,„ i . , , j r°° du (m)(u) / . , , 

/ — ^ A ~ L du (I - u)^ 1 sm(ujtu )=t~ 1 / — ^ i^J.[i-S al coB(iJt) + 

J0 7T U) Jo J TT U) 1 V 

-(a - 1) f 1 du (1 - n)- 2 cos(u,tu)) ~ r 1 f° ^1^0 (52) 

where we have neglected in the last expression the oscillatory integrals by invoking t S> Tq. By 
inspection of (|C8p we conclude that for m > 3 the frequency integral in (f23j) can be approximated 
byt^l/^drr^r)!. 

For an Ohmic bath m = 1, the last integral in (|52p diverges. We now argue that one can replace 
(JRh)(w) by 0i)(O) = 27 (/i/?)" 1 on the left-hand side of (JS2J, which becomes 

/ — i —I du (1 - u) - 1 sm(ujtu) = ± <±-L (53) 

Jo -k 10 Jo 2a 

in the limit t ^> Tq. (We have used J duJsin(ujtu)/LO = tt for tu > 0.) Note that this amounts 
to replacing $lh(t) by a white-noise correlator 2^/(hf3)~ 1 5(t) in (|42p . Let us estimate the error 
introduced in the frequency integral in (|43|) by this substitution. This error is given by the left-hand 
side of ([52]) modulo the replacement of (9ft/i)(a>) by (3Rh)(u)) — (5ft/t)(0). Disregarding oscillatory 
integrals as in the case m > 3, the error is equal in the limit t ^> Tg to t^ 1 J °° do; ((K/i)(w) — 
(K/i)(0))tJ -2 /7T. The latter integral converges since (3ft/i)(u>) — (3ft/i)(0) behaves like to 2 for small u>. 
Comparing with ([53|) (see also (|C9p ). one concludes that the relative error introduced in (|4"3"]1 by 
the substitution of (3ft/i)(w) by its value for to = is small, of the order of Tg/t. Hence, for m = 1 
the frequency integral in (|4*3"P can be approximated by (3ft/i)(0)/(2a) = a~ l J °° dr 5R/i(r). 

Collecting the above results and integrating (|4^|) with respect to time, we find in the Ohmic 
case m = 1 

Dr k (s, s ') = ' 

''decvS, S ) — 

and in the super-Ohmic case m > 3 

D t pcak (s,s') = (- — J—y:) (super-Ohmic) , (56) 

, ( ,, _ ( 2{B 2 )h 2 (i 2 \^ t ent (s,s>) 

t dec {s,s) - { Ca{SjSf) 2 lJo °° dTTmT) \) rjL/a ^ 

with the proviso Tg <C £dec( s ) s ') ^ Ts,Tp. We can interpret the growth of D pea like t 2a+l in the 
Ohmic case by saying that for fixed x and x', in the Markov regime Dt must be proportional to 
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(Ohmic) , 


(54) 


(2a + 1) (£ 2 ) ft/? \a-+i /i e nt(s,s')y^ 


(55) 



t(x la — x a ) 2 (the fact that D t oc t is well known [3J]); the indicated time behavior of Df e then 
follows by replacing (x, x') by (ets,ets'). 

By using \h(r)\ < (B 2 ) and Rh^r) ~ for r > T B (Sec. IIIIAj) . one finds that the integrals 
Jp 00 dr 3ft/i(r) and | J °° drr 3ft/i(r)| are at most of the order of (B 2 )Tq and (B 2 )Tg, respectively. 
If \s'/s — 1| is not close to unity (so that c a (s,s') in (|47h is not very large) the factor inside the 
parenthesis in (|57p is of the order of (hfi/Tg) 2 or larger. Thus, for coupling strength r\ < h(3/Tg the 
condition t C nt(s,s') < idec^s') holds in the Markov regime for super-Ohmic baths. The situation 
is different for Ohmic baths: then, by (|55p . the condition in question is violated even for small rj if 
the entanglement time t en t(s,s') is large enough compared with h/3. More precisely, still assuming 
that c a (s, s') is of the order of unity, £dec( s > s ') becomes smaller than t en t(s, s') when 



t e nt(s,s') . {B 2 )k[i 



> 



(58) 



h/3 ~ ri 2 f™drm(T) ' 
For super-Ohmic baths, the decoherence time ()57[) decreases by increasing a for 77 < hp/Tg 
and |s'/s — 1| not close to unity, i.e., provided that td ec (s,s f ) > t ent (s,s'). Then t^ l^dec ~ 
(??Tg/(/i/3)) 1_1 ' a < 1. Thus, for fixed weak enough coupling strength rj, nonlinear pointer-bath 
couplings always win over a linear coupling in efficiency for decoherence. This is in striking contrast 
with what happens in the Ohmic case. Actually, for a Ohmic bath nonlinear couplings become 
less efficient than a linear coupling when t en t(s,s') is large enough so as to fulfil (|58p . More 
precisely, we find by using c a (s,s') as 1 and (f55j) that the decoherence time ([55]) is larger in the 
nonlinear case a > 1 than in the linear case a = 1 by a factor r£ Adec" °^ ^ ne order of 

Finally, it is worth mentioning that Ohmic baths win in efficiency over super-Ohmic baths. This 
can be shown by noting that (idee /^dec m ) 2a+1 is equal (up to a numerical factor of the order of 
unity) to the product of | J °° drr Rh(r)\(h/3 J °° dr ^^(r))" 1 by H(3/t^^ m . Since the last factor 
must be small compared with 1 for consistency (recall that hf3 < Tg), it follows that ^^(s, s f ) is 
smaller than i^ m (s,s'). 

One may wonder if the results of this section could be strongly modified if a direct coupling 
between the object S and bath B (which we do not admit in the present model) was allowed for. 
It is clear that one can answer this question by the negative when the object-pointer coupling 
constant e is large enough, i.e., for small enough t e nt- I n order to estimate how small must be 
tent) let us couple S and B via the Hamiltonian H$b = A a (S/5s) a B. This Hamiltonian has 
a magnitude comparable with the pointer-bath coupling i|3]) in the initial state (|4a[) . We first 
consider Ohmic baths. It is known that the decay of the off-diagonal matrix elements (s\ps(t)\s') 
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resulting from the coupling Hge then goes like exp(— £/Td ec (s, s')) in the Markov regime (we ignore 
here the object-pointer coupling) [3]. If T^ ec (s,s') <C T$, a condition fulfilled e.g. if [S, Hs] = 
(pure dephasing regime, T$ = oo), the corresponding decoherence time is given by T^ cc (s,s') = 
h 2 (5s/A) 2a \s' a — s a \~ 2 / f dr 3ft/i(r) [2l|]. The ratio between Td ec (s, s') and the decoherence time 
d5SD for s' = s + 5s is [{c a 7])-H-^{h(3) 2 {B 2 )/ / °° dr3fJ/i(r)] 2Q! /( 2Qf+1 ) up to an irrelevant factor. 
Taking into account that J^° dr K/i(r) < (B 2 )Tg, we see that it is well justified to neglect the 
coupling of the object with all degrees of freedom of the apparatus but the pointer provided that 
tent / {h(3) <C (h(3/Tj3)(c a rj)~ 2 . For a super-Ohmic bath, if [S,Hg] = then the modulus of the 
off-diagonal matrix element (s\ps(t)\s') decays to a nonzero value under the coupling Hsb (for a 
discussion on this saturation of decoherence see e.g. [48fl). whereas indirect decoherence via the 
pointer leads to a complete decay of the object-pointer coherences (this decay being given by 
the decoherence exponent ([56]) ). It is also easy to show that T^ ec (s,s') is much larger than the 
decoherence time flEE} provided that t cnt /(h(3) < (H/3/T B ) 1+1/a (c a ri)~ 2 ~ 1 / a . 

D. Bath of harmonic oscillators linearly coupled to V 

To study the transition between the limiting time regimes discussed in the two preceding sub- 
sections, let us consider a bath of N S> 1 harmonic oscillators, Hq = Yl,u ^>u(blb u + 1/2), coupled 
to the pointer via a coupling agent B linear in each of its creation and annihilation operators b v 
and b u , B = ^2 u {K u bi + K u b v )/y/N [24|]. Here uj v is the frequency and n v the coupling constant of 
the vXh. oscillator. We shall take the following specific choice for the power spectrum function: 

N 

J(u) = ^ E K\ 26 ^ - ^) = 7^ m e-" 2/ ^ (59) 

wherein m is an odd positive integer, 7 > 0, and Up is a cutoff frequency. We recall that the case 
m = 1 corresponds to an Ohmic damping, whereas one speaks of super-Ohmic damping for m > 1. 
For instance, m = d or d + 2 for a phonon bath in d dimensions, depending on the underlying 



symmetries 25J. As is well known 2J, |25J, the imaginary part of the bath correlation function 
h(t) is temperature-independent, its Fourier transform being given by i(Qh)(uj) = J(u>) for uj > 0. 
By the KMS property (JCH) this implies 0i)(u) = coth(Hf3u/2)J(\u\). If w D = Huj d (3 > 1, the 
thermal time Tg = fh/3 is the largest decay time of ?R.h(t). The other time scale characterizing the 
variations of $lh(t) is the inverse cut-off frequency tg = uj^ < T$. By (jM]), the decoherence and 
entanglement times in units of Tg, T& ec = t& ec /TQ and r ent = t en t/Ts, are given by 
^ent _ Jo°° dw coth(tu/2) w m e~ w2 l w l | J Tdcc dr - 



c 2 V 2 2 J °° dw coth(w/2) w m e- w2 / w l 
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(60) 



where we have expressed $lh(t) in terms of its Fourier transform and relied on (|59p . We did not 
write explicitly in ([UUjl the dependence of T ent , Td ec , and c a on (s, s') . The right-hand side of (fB0|) 
is shown in the inset in Fig. [3] We have computed numerically the integrals appearing in this 
right-hand side for various values of a, to, and wd, so as to obtain Td cc as a function of r en t and r/. 
The main results are shown in Figs. 0and|31 For fixed a and rj, the plain curves representing Td ec in 
Fig. [3] split by increasing r cnt into distinct branches corresponding to distinct m's, as predicted by 
([33]) and ([36]) . This splitting occurs when r^ ec is in the transition region w^ < Td cc < 1. After this 
splitting Tdcc is larger for larger m. In particular, a Ohmic bath (to = 1) has a smaller decoherence 
time than a super-Ohmic bath (to = 3,5...) as stated above. For comparison, the power law 
behaviors found in Sec. IVIBI and IVICI in the small time (rd ec <C w^ ) and Markov {t& cc 3> 1) 
regimes are also shown in Fig. [3] (broken lines). A remarkably good agreement between the exact 
and asymptotic behaviors of t^ cc is obtained: the exact results are well approximated by their 
small-time behaviors (|5ip up to Td cc < w^ and they are hardly distinguishable from the Markov 
approximation as soon as Td ec > 1- Our aforementioned statement that a nonlinear pointer-bath 
coupling is more efficient for decoherence than a linear one when r en t is not too large (and even 
for arbitrarily large r en t if to > 3 and rj is small enough) is well confirmed. Indeed, it is seen in 
Fig. H] that for a pointer-bath coupling strength rj <C 1, t^ cc becomes significantly smaller when 
the value of a is increased from a = 1 to a = 3. If the dotted lines (Markovian results) in Fig. [3] 
were drawn farther to the right, the two lines corresponding to (a, to) = (1, 1) and (a, to) = (2, 1) 
would intersect; after this intersection (not shown in the figure) the reverse situation of higher 
values of a leading to higher values of Td ec occurs. In contrast, for to = 3, 5, . . ., the dotted lines 
associated with a = 1 and a = 2 never intersect (their are parallel); hence Td ec decreases with a 
and r en t < T dcc for all values of r ent (more precisely, r ent < lCT 4 T dcc for a = 1 and r ent < 10 _2 r dec 
for a = 2). We also emphasize that t^ cc increases in Fig. [JJ with the cut-off frequency ojq- Even 
though the results in Figs. [3] and H] correspond to the simplifying choice of a bath of harmonic 
oscillators with power spectrum function (|59p . for more general baths they should still give the 
correct qualitative picture. 



E. Bath at very low temperature 

We have so far considered baths at finite temperature. Motivated by experiments in solid 
state physics, we shall now discuss the case of a bath initially in thermal equilibrium at very 
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a = 2 



a = 3 ^L--" - - 



1.5 L 



FIG. 4: Decoherence time t^ cc in units of Xg = H/3 as a function of the pointer-bath coupling strength r\ 
for the same bath as in Fig. [3] with r en t = 0.1, m = 3, and c = 1. Three distinct values of iba are shown: 
iu,d = 2 (solid curves), wd = 5 (dotted-dashed curves), and iud = 10 (dashed curves). For each of these 
values, Tdcc is shown for a = 1, 2, and 3 (from top to bottom). 

low temperature. Strictly speaking, for the equilibrium apparatus initial state (J4bj) extremely 
low temperatures have to be proscribed because of our hypothesis h/3 <C Tp. However, taking 
e.g. T-p = 1 s, this separation of time scales holds even for the smallest temperatures that can 
be achieved in experiments. Furthermore, the stability conditions (|40p have a better chance to 
be met at low temperature T since A t h decreases with T. To be specific, we consider the same 
bath of harmonic oscillators as in the previous subsection, but now in the limit wo = tkonP 3> 1. 
Then only spontaneous emission plays a role in the pointer-bath interaction. In other words, 
iyRh){uo) = coth(hf3uj/2)J(\uj\) can be approximated by J(|u;|). The zero-temperature variance of 
B equals (B 2 ) = J du; J(\uj\)/(2it). For our choice (|59p of the power spectrum function, this gives 
2tt(B 2 ) = juj™ +1 ((m - l)/2)!. The analog of flBDJ reads 



(co D t 



ent; 



\2a 



1 



cIVd 



<^D< dc 



, , dvv m e~ v / duu a e~ ivu (61) 

((m-l)/2)!J Jo V ' 

where tjd is now the pointer-bath coupling strength in units of Hu>d, t]d = (B) 1 ' 2 A a /(hu>D)- 

Eq. (f61~j) holds provided that uj^tdec *C wd, i.e., td ec ^ fifi. This equation is the same as ([60]) apart 

from the substitutions T ent — > uJDt en t, Td ec — ► WDtd ec , V ~> VD and coth(w/2)e~( w > WD ' > — ► e~ v . 

Explicit formulae for td ec can be given as before when td ec is small or large compared with ui~j^ . 

One reads the small-time result directly on (|5ip by transforming this expression according to the 

recipe mentioned above. This gives u)otdec oc (c Q T) D u}Dt e nt) a for tdcc *C u>p . Similarly, 

for a super-Ohmic bath tdec oc c a Vd *ent when uj^ <C td oc <C /i/3. To find the proportionality 

factor, it is enough to realize that | / °° d.TT$th(r)\ = J °° dw (3f?/i)(cj)o'~ 2 /7r has to be interpreted 
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W D t 6 



FIG. 5: Decoherence against entanglement times in units of lh^ 1 for the bath at zero temperature of Sec. lVIEl 
with a = 1, r\B = 1, and m — 5, 3, 1 (solid curves, from top to bottom). The approximations (|5ip and (|63[) 
for wijtdec <C 1 and Wuidec 3> 1 are shown as dashed lines. 

in (|57p as 7W^) _ / °° dv v m ~ 2 e~ v /it = 2w^ (B 2 )/(m — 1). The Ohmic case requires some extra 
work. For indeed, replacing (JR.h)(uj) by J(|w|) leads to a vanishing second member in (|53j) . even 
though the frequency integrals in the second and third members in (|52j) are still divergent. One 
actually finds 



OO pi 

2 



due * / du{l — u 

ii Jo 



) sm(u>Dtuv) ~ as u>ot — > oo (62) 

WD* 



with k\ ~ 0.2886 and /c Q = fei — 1 — 1/2 — ... — l/(o — 1) if a > 2. Substituting the frequency 
integral in (|43p by the right-hand side of (|62p . one gets 



(63) 



icnt(s,s') - (c Q (s,s')r? D ) t dec (s,s')^ln(u; D t dcc (s,s / )) +^a - ( 2 «) ^^ (Ohmic) 

1/q 



\ I/O 

c a (s,s')r] D /Vin- 1) t dec (s,s') 



(super-Ohmic). 



Instead of going through a proof of (j62|) . which would lead us too far into technical details, let us 
compare the formulae (|63H to the exact results obtained by numerical evaluations of the integrals 
in (I6ip . It is seen in Fig. [5] that the approximate values (|63h closely follow the exact curves when 
uj£>tdec becomes large (in fact, even for uotdcc — 2 in the Ohmic case m = 1). Similar pictures 
are found for higher a's. Let us remark on (|63p that for a given t ent , the ratio between the 
decoherence times for Ohmic and super-Ohmic baths is logarithmically small in the dimensionless 
time uJotdcc- Hence a Ohmic bath is not dramatically more efficient than a super-Ohmic bath at 
very low temperature, in contrast with our previous findings at "high" temperatures. 
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VII. DERIVATION OF THE RESULT 

We here fill in the derivation of the results presented in Sec. [V] An alternative derivation of 
([32j) based on the time-dependent Redfield equation can be found in [49]. Our approach below is 
non-perturbative in the pointer-bath coupling but makes use of the QCLT (Sec. IIII A[) which holds 
due to the additivity (|3|) of the bath coupling agent. 

A. Object-pointer dynamics 

According to the results of Sec. IIIII we may drop the Hamiltonians H$ and H-p in the full 
Hamiltonian H, with the proviso that the object initial state ps is replaced by p s {t) given by (|16p . 
This means that at times t <C T$,Tp the exact evolution operator of S + V + B in (JSJ) can be 
approximated as e ' itH / h ~ e - [tHB ^ h W(t, 0) e -'<H s +H v )/h^ w h ere 

W(t,0) = jtH B /h e -it{H B +H vs +H VB )/h = Texp |_ i f dr ( eSP + x a B(r))\ (64) 

is the approximate evolution operator in the interaction picture. Note that B(t) is different from 
B as soon as r > tg. Since we do not assume here t to be small compared with ts we keep the time 
dependence of the bath coupling agent in (|64l) . In view of the product structure p(0) = ps <8> Pvb 
of the initial state and by cyclic invariance of the trace, the object-pointer state (JSj) becomes 



Pvs (t)~tT B (W(t,0)p s (t)®p VB W(t,0)^ , t<^T s ,T v . (65) 

The pointer Hamiltonian H-p is absent in (|65p since e ~ ltHv ' h p-pjge ltHv ' h ~ p-p& at times t <C T-p 
for the initial states under study, see Sec. IIIB 

The approximate evolution operator (|64h can be simplified by using the exact identity 

W(t,0) =e-' lteSP/h Texpl-^ f dr(X + reS) a B(r)\ . (66) 

We forego the proof of this (generalized Baker-Campbell-Haussdorff ) identity, which uses the role 
of the momentum as generator of displacements, e -lteS ' p ' fi '(X + teS) a e i ' h = X a . Employing 
([55]) in ([US]) and setting x s (t) = x — tes and x' s , (t) = x' — tes' as before we get 

( S ,x\pp S (t)\ S ',x') = (s\p° s (t)\ S ')( Xs (t)\tvJrexp!.-^J dTX s (t-T) a B(T) 



Tex V {-^f\rx' s ,{t-rrB{T)\ W,(t)) . 



(67) 

The next step consists in evaluating the trace over the bath in (j67p by taking advantage of Wick's 
theorem ([T2]) . We discuss the two initial states (|3a|) and (j4b>[) separately. 
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B. Partial-equilibrium initial state 



For the partial-equilibrium initial state (|4a|) one has p-pg = p-p p^' and the last matrix 
element in (|67p is the product of a pointer and a bath expectation values, 

(s,x\ Pvs {t)\s\x') = (s\p° s (t)\s') {xMpvWAt)) K t{x S (t)Ai(tY,s,s') (68) 

with 



K t (x,x';s,s') 



Texp 



h 



Texp 



dTx' s ,(-T) a B(T) 



dTX s (-T) a B(T] 



t 



(69) 



Here (•) = Zg tr^-e P"B) denotes the average with respect to the free bath thermal state. The 
QCLT and the additivity (J3]) of the bath coupling agent imply 



Ft,o[k,t\ 



Texp<J— - / dTk(r)B(r) 



t 



Texp<J-- / dTl(r)B(r) 



<X!) i "^2 J dTl J' d M k (n) - l(Tl))(k(T 2 )h(T 2 , Tl ) - l(T 2 )h(n,T 2 ))\ , (70) 

where k(r) and 1{t) are two arbitrary real-valued functions and h(ri,T 2 ) = h{r\ — t 2 ) is the two- 
point bath correlator, see (jl(J|) . The identity (|7Up is equivalent to Wick's theorem (|12p . For ordered 
times t > t\ > t 2 > ■ ■ ■ > t n one actually gets (fT2j) from (|70|) by setting k = in ([TO]) and taking 
the functional derivative of both members with respect to l(ti), . . . , l(t n ) at / = 0. The proof of 
converse statement is deferred to Appendix [Dj By using the parity properties ?R.h(T) = 3ft/i(— r) 
and 9/i(r) = —Qh(—T) of the real and imaginary parts of h and employing (|70p in (|69p we get 



K t (x,x'; s, s') = e~ 



-Dt(x,x';s,s')— i</>t(x,a;';s,s') 



(71) 



with a decoherence exponent L>j and a phase 0t given by 

1 



{D t + icj> t ){x,x';s,s') = ¥ J dn / dr 2 ( 3 ;' s ,(-T 1 ) a - ^(-n) )^ (x' fl ,(-r 2 ) a - x s (-T 2 ) a ) 



xm( n - r 2 ) - \{x' s ,(-t 2 y + x s (-T 2 rph( Tl -T 2 )\. 

Thus (|68|) reduces to the result ([32]) announced in Sec. IV Al 



C. Equilibrium apparatus initial state 

Before deriving the expression corresponding to (f68j) in the case p-pg = p^g , we determine the 
pertinent initial density matrix of the pointer, 

R (x,x') = (x\tv B (p { ^)\x'). (72) 
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It is convenient to introduce the x-dependent bath average 

(0 B ) x = Z-Hv B (0 B e-^ H ^ a ^) , Z a . = tr B (e-^+^)) (73) 

with x a real number (the pointer position, for us). Note that 

e -l>(H B +x"B) =e -l3H BTexp t^ r dz B(-iz)\ (74) 

with B{— iz) = e z B ' Be~ z B ' . The normalization factor Z x can be determined by applying 
Wick's identity ([70]) with t = -ifr/3, k(r) = 0, and /(r) = x a . This gives 

f x 2a f3 l0 ] 
Z x = Zq exp <^ — - — \ (75) 



To = T3 / dz i / dz 2 ^(-^2) • (76) 

bP Jo Jo 



with 

ftP Jo Jo 

By using the analyticity and KMS properties of the bath correlator h(r), one can show that 
< 70 < fr(3(B 2 )/2 and that 70 = 7(0) coincides with the following integral evaluated at t = 

7 (t)= J dr3/i(r) . (77) 

J —00 

Details of this derivation are deferred to Appendix [Cl 

Replacing (|75p into the approximation (|17p for the apparatus initial state and inserting the 
high-temperature expression of (x\e~^ Hv ' 2 \y) (see Sec. [TTJ) yields 

flops') = ^ J dyZ y (x\e-^ H ^ 2 \y)(y\^ H ^ 2 \a/) 

= ZoZ-pb' 1 /'dye^ 7oy2Q/fi e"^ (y(x)+y(x ' )+2y(y))/ S" 4 ^ 2((x ~ y)2+(x '^ )2)/A ?h (78) 

with A t h = 2irh((3/M) 1 ' 2 . The stability condition (|40p does not guarantee that V(x) compensates 
—jox 2a /h when x — ► ±00. If this is not the case, i.e., if the effective potential V e g(x) = V(x) — 
70 x 2a /h has the shape shown in Fig. [21 the integrals in (|78p diverge. This reflects the fact that the 
pointer interacting with the bath will tunnel to infinity after a certain time. Since we restrict our 
attention to initial states describing a pointer initially localized inside the potential well of V e s(x), 
we shall disregard this convergence problem by adding to V(x) a positive potential vanishing 
for \x\ < W e ff and diverging exponentially for x — ► ±00. This regularization trick amounts to 
replace p^g in (|72p by the local thermal state of the apparatus discussed in Sec. IV Bl After 
this regularization, the main contribution in the y-integral in (|78p comes from small values of 
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y, \y\ < Ath- In fact, the first exponential in this integral is a slowly varying function on the 
scale A th since (/37o//i)~ 1/(2a) > 2 1 ' a A th > A th , as follows from flU and 7o < h(3(B 2 )/2. Due 
to the presence of the last exponential in (|78p . this first exponential can be approximated by 
e /37o(x+z ) a /(2 a h) an( ^ t a k en ou i f the integral. Similarly, the second exponential varies noticeably 
on the scale A t h 3> A t h and can be approximated by Q-P( v< y x ) +V ( x )+zv(x/2+x /2))/4 anc j t^en ou t 
of integral in ([78]). Thus (gOD and A th > A th entail 

Ro(x, x') ~ Z^s e -/3(^ ff (-)+K ff (-'))/2 e - 27r2 (.-,')VA t \ (79) 

with V e s(x) given by (|35p and Z-p^g = f dxe~^ eB ^ x ' . We have used in (|79p the approximations 
/3F(x/2 + x'/2) ~ /9(F(s) + V(x'))/2 and /^(x + x') 2a /(2 2a H) ~ /37 (x 2q + x' 2a )/(2/i). This 
introduces an error which is negligible against (x — x') 2 /A 2 h for |x|,|a/| < Ath an d A t h 3> Ath- 
Hence the pointer is in a Gibbs-type state with an effective potential V c s(x), as announced in (|38p . 
We can now proceed to evaluating ([67]) . Repeating the steps yielding to ([78]) and using the 
notation ([75]) . 

< S ,x|pp 5 (t)| S ',x') = Z"i <fl|^(t)| S ') /"dyZ^ar.Wle-^/^Xyle-^/VyC*)) x 



t r s /■* 



Texp j— i / drx',,(t - r) a 5(r)| Texp |— i / drx s (t - r) a 5(r)| 



(80) 



a 



We set 6B(r,y) = B{t) — (B{r)) y and consider the (quantum) characteristic functional 

t f , /-t 



Texpj-^y dTfc(r)<5S(r,y)} TexJ-^J drl(r)SB(r,y)\ 



(81) 



;</ 



It is shown in Appendix [D] that all the correlation functions (SB(ji,y) ■ ■ ■ 5B{r n ,y))B,y are inde- 
pendent of y, i.e., 

Ft,,, [M] =**,<> [M] (82) 

for any y, t, fc(r), and l(r). Wick's theorem (fi~2]) also entails (see Appendix ID"]) 

~ 2r Q 

(B(r)) x = -— 7 (r) (83) 

with 7(r) given by ([77p . Formula (|83[) is reminiscent of linear response theory since j(t) = 70 — 
h J dr x(t) @(t — T)/2, where 0(r) is the Heaviside function and x( T ) = ~(2/K) 9(t) Qh(r) the linear 
susceptibility. Let us point out, however, that ([83]) is exact to all orders in x a . Collecting the above 
results one finds 

(s,x\pvs(t)\s',x') = (s\p° s (t)\s') R t (x s (t),x' s ,(ty,s,s') Kt{x 3 (t),x' s ,(t);s,s') (84) 
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where K t (x, x'; s, s') is given by ([69]) . 

R t (x,x'; Sl s') = Ce-^(*) +v *We-^^)^x 

J d£ ex V {-e + V^Z^^-2ie 9t(x,x'-,s,s')} , (85) 

C is a time-independent normalization constant, and 

9t(x,x';s,s') = (8vr 2 )-f^ j* ' dr 7 (r) (x' s ,(-r) Q - x s (-r) Q ) . (86) 

Therefore, dTTJ) and (jEU) account for d37j). Moreover, by using (Jib]), ([16]), ((TJj), (JHH), and X = 1 
one easily establishes that Ro(x, x'\ s, s') = Ro(x, x'). Taking t = in ([85]) . evaluating the Gaussian 
integral and comparing with (j79j) . one gets C = n~^Zp eS . 

D. Justification of the approximation (|38p for i?o 

We here want to derive the inequality 

gt(x, x'; s, s'f < (87r 2 )~ a ^P^ D t (x, x'; s, s') . (87) 

Let us first point out that, putting together (|87p . the stability condition (|40|) . the separation of 
length scales A t h <C A t h, and the bound 2{3yo/h < rj 2 h /A 2 £, ^ follows that gt(x,x'; s,s') 2 <C 
Dt(x,x';s,s') uniformly for all (x,x') and (s,s'). This explains why the general expression (|85f) 
reduces to <[38|> for short times t satisfying Dt(x,x',s,s') < 1; then \gt(x, x'; s, s')\ <C 1 and the 
phase factor inside the integral in (|85p can be neglected; by performing the resulting Gaussian 
integral, one gets Rt(x,x'; s,s') ~ Rq(s,s'). In the special case a = 1, the integral in ([85]) can be 
evaluated exactly for all times i. This leads to Rt = i?o e _1< Pt. Here and in what follows (f>' t , 

(f>1, etc, denote real phases irrelevant for decoherence. Replacing the latter value of Rt into ([57]) . 
the factor e~^ 9t ' can be dropped by invoking gr| <C Z^ again. We are thus led to 

e- Dt{x ' x '' s ' s,) R t (x, x'; s, s') ~ e^*'*'"' - '^'^^, z') (a = 1) (88) 



which is now valid for all times t <C Zs , T-p . The integral (j85p can be evaluated exactly for 
a = 2 as well. By (j87[) and the same arguments as above, for a > 1 the stronger condition 
gt(x,x';s,s') 2 <C (A t h/A t h) Dt(x,x';s,s') holds. Using also the restriction \x\, \x'\ < A e g « A t h 
coming from the factor in front of the integral in (1851) . one obtains for all times t <C Ts, T-p 



p -D t (x,x';s,s')-i(j>' t " r> I A 

e-W^Rtix, x>; s, s>) x \ , ° ^ ' (a = 2). (89) 
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Notice that this equation is consistent with Rt — Ro at times t satisfying Dt <1. 
Proceeding towards the inequality ([57]) we rewrite ([55]) as 

A 2a gt(s,s;s,g) 2 = (/ ^^^J drcos(wr)(x^(-r) a -x s (-r) a ) 



where 7(0;) = 7(— w) > is the Fourier transform of 7(i), see ([771) . By using 70 = 7(0) = 
f du)'j(u>)/2ir and the Cauchy-Schwarz inequality, 

(Sn2) ^gt(x,x';s,s'f < 7o|°° ^l{oj) (^dr cos (u,r) (^, (-r)" - z s (-t) q )) (91) 



\ 2 h 



The integral over u in the right-hand side of (|9ip can be bounded with the help of (|C11|) by 

(92) 



Thj r ' U ' "'"-' 1 [«/ ,t dre- iwT (x / a /(-r) a - a j 8 (-r) a ) 



2 ./_oo 2?r 



Comparing (|92[) with (|34f) . we bound the last quantity by h?f3Dt(x,x',s,s') and have thus estab- 
lished the inequality (f87|) . 



VIII. CONCLUSION 

Let us summarize the main results of this paper. We have investigated a model for a quan- 
tum measurement in which the entanglement produced by the interaction between the measured 
quantum object and the pointer is simultaneous with decoherence of distinct pointer readouts; 
the apparatus (pointer and bath) is taken initially in a metastable local thermal equilibrium, not 
correlated to the object. Our model has four parameters: the object-pointer coupling constant 
e, the thermal variance {B 2 ) of the bath coupling agent, the temperature T = (/c^/?) -1 of the 
bath, and the exponent a in the pointer-bath Hamiltonian ([3]). One may construct out of the 
first three parameters two relevant dimensionless constants. The first one is the entanglement time 
Tent = A(e 5s fr(3)~ l in units of the thermal time h(5. Here 5s is the separation between neighboring 
eigenvalues of the measured observable and A the uncertainty in the initial pointer position. That 
entanglement time r en t describes the efficiency of the pointer-bath interaction (a coupling is effi- 
cient if Tent is small). More precisely, r en t is the time after which pointer positions corresponding to 
distinct eigenvalues s begin to be resolved. The second dimensionless combination is the coupling 
energy r\ = (B 2 ) 1 ' 2 A a f3 in units of ksT, which measures the strength of the pointer-bath coupling. 
We have found that, after a certain time id ec , the object-pointer state is close to a statistical mix- 
ture of separable states ^ s p s \s)(s\ (8) pj,, with p s = (s\ps\s), ps the object initial state, and pf, a 
distinguished pointer state depending on s. The decoherence time idee = ^^dec <S T$,T-p is given 
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by (we ignore here numerical factors, given explicitly in (f5Tj) . ([55]) . and (|57|) ) 



'',[«-<■ -v I?? 1/a ^nt) , 7 



-£ri if ^dec ^ *B (interaction-dominated regime) 
2a if tdec > T B for an Ohmic bath (Markov) (93) 



2cM-l 

1 if *dec ^ ^jg for a super-Ohmic bath (Markov). 

For reasonably strong pointer-bath coupling and not too strong object-pointer coupling, the deco- 
herence time idee (needed for the transformation of linear superpositions into statistical mixtures) 
can be so small that the whole measurement is performed without producing a Schrodinger cat state 
as an intermediate step. Two distinct regimes ought to be identified in (|93|) : in the interaction- 
dominated regime, idee is shorter than the characteristic time ig after which the bath correlation 
function h(t) differs significantly from its value (B 2 ) at t = 0; in the opposite Markov regime, 
one must wait more than the bath correlation time Tg, i.e., the largest decay time of h(t), to 
obtain the required statistical mixture. While idee presents a universal behavior in the interaction- 
dominated regime (it depends on the bath through the single parameter rj), in the Markov regime 
it is determined by the small- frequency behavior of Qh(t), (Qh)(ui) ~ —i^yuj m . Larger values of 
tdec ar e found for larger m's, with a significant change of behavior between m = 1 (Ohmic bath) 
and m > 1 (super-Ohmic bath), see ([93]) . In both time regimes, idee strongly depends on the 
nonlinearity exponent a, as illustrated in Figs. [3]and|31 Smaller decoherence times are obtained 
for larger a's save for the Markov regime if m > 1 and r\ > hfi/Tg or if m = 1 and i] 2 T CIlt > hft/Tg, 
where the reverse statement holds. The linearization of the pointer-bath interaction with respect 
to the pointer position (dipole approximation) may then lead to an over-estimation of tdec in the 
interaction-dominated regime or for super-Ohmic baths in the Markov regime. For a bath at very 
low temperature, (f93j) still holds with Td ec and n replaced by tdec/^B & n d r\D = (B 2 ) 1 ' 2 A a tQ/h, 
save for the Ohmic case where idecAent becomes logarithmically small in tdec/^B- 

Several generalizations of our results may be of interest. The first one concerns measurements 
of observables with continuous or dense spectra. One must then allow for a finite resolution Ss in 
the measurement result. Unlike in the case of discrete non-degenerate spectra studied in this work, 
the "final" object-pointer state will not be a separable state because coherences for pairs (s,s') of 
close eigenvalues (\s — s'\ < Ss) are damped on a smaller time scale exceeding the time duration 
of the measurement. A second generalization concerns the bath, assumed in this paper to consist 
of independent degrees of freedom. As stated in Sec. HH it can be shown that the validity of the 
QCLT extends to baths of interacting degrees of freedom if the correlator {B^B P ) decays more 
rapidly than l/\fj- — v\ (see 23J for a related version of the QCLT in this context). This implies 
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that our results apply to a broad class of baths including certain interacting spin chains. 

It would be interesting to investigate concrete models for the object and pointer involving 
projective measurements in the "no-cat" regime (decoherence fast compared with entanglement), in 
connection with recent experiments in solid state physics. We should also mention that our results 
can be of interest in a broader context. Actually, we have studied quantitatively a new scenario for 
decoherence. In this "indirect decoherence" scheme, the decay of the quantum coherences of the 
small system (object) does not result from a direct coupling to the many degrees of freedom of a 
bath, but rather from a strong coupling to few degrees of freedom of the environment only (here, 
to the pointer). These few degrees of freedom are in turn coupled to all others bath coordinates 
and serve as an intermediate in the decoherence process. 
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APPENDIX A: APPARATUS AS A NONINTERACTING INFINITE GAS 

Let us consider a gas of N = N + 1 noninter acting particles with mass m v , momentum P v , and 
position X v (y = 0, . . . , N) submitted to a slowly varying external potential V{x). To simplify the 
discussion we restrict ourselves to a one-dimensional geometry. Let M = ^2 U m u and P = ^2 u Pu 
be the total mass and momentum and X = ^2 V m v X v jM the center-of-mass position; R u q = 
A/" -1 ' 2 (X v — Xq) and P u q are the relative positions and their conjugate momenta. Expanding 
V(X„) as V{X) + {X V -X)V'(X) and using A^-X = V^^o-E^AoW. the Hamiltonian 
of the gas reads 

P 2 N / Kf \ N P 

H ™ = h + Nv{x) + v^£ (i - ^ W w + £ sfefl-fl. (ad 

ffp S V ' v ' 

where •£ is the N x N matrix with inverse ^ _1 = (5^ u — m v /M)^ v=l . The pointer V is the center- 
of-mass degree of freedom. Its Hamiltonian H-p is given by the two first terms in (|A1|) . The 
bath B is constituted by the N relative degrees of freedom. Its Hamiltonian Hq is the last term 
in (|Aip . The third term in (jAip describing the coupling between V and ,8 has the form ([3]) if 
V(sc) = N~ x (a + l) _1 x a+1 and B is given by © with B v = (1 - Afm v /M)R u0 . If the measured 
system is strongly coupled to the total momentum P of the gas, one obtains a tripartite model of 
the kind discussed in Sec. [H] (although Hq does not satisfies all our hypothesis). 



APPENDIX B: APPROXIMATION FOR THE APPARATUS EQUILIBRIUM STATE 

In this appendix we justify the approximation (117p for the initial Gibbs state p^Q of the 
apparatus. Moreover, we show that Pp(t) = e~ ltHr ' h p-pe ltHv ' !i ~ pp when t <C T-p for the quasi- 
classical pointer states p-p considered in Sec. HH A similar result holds for p^pg • 

We recall that Tp is defined by Tp = (M/V"(0)) 1/2 . Taking V(x) ~ V"(0)x 2 /2 and invoking 
©, ©, and (jSJ one easily finds that tip(X(t) 2 pp) - A 2 and trp(P(i) 2 p-p) - Ap 2 are equal to lowest 
order in time to (-A 2 T p 2 + Ap 2 M- 2 )t 2 « A 2 T v 2 t 2 and (-Ap 2 T v 2 + y"(0) 2 A 2 )i 2 « Ap 2 T v 2 t 2 , 
respectively. Hence Tp can be identified with the time scale for significant evolution of X(t) and 
P(t) when the pointer is in the quasi-classical state ([7]) (Sec. IIII A[) . One easily convinces oneself 
that pp = pp(X,P) is an operator- valued function of the position and momentum operators X 
and P. Letting pp evolve under the Hamiltonian H-p up to time t amounts to substituting X 
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by X(t) and P by P(t), see @. This shows that p° v (t) = p(X(t),P(t)) ~ p? as X(t) ~ X and 
P(t) ~ P for t < Tp. 

In order to approximate p^g we shall rely on the Baker-Campbell-Haussdorff formula 

e A Q C = e A+C+[A,C]/2+[A,[A,C]]/12+[C,{C,A]}/l2+- ^y 

wherein A and C are any two operators. After a few transformations, (|Bip becomes 

e W = e 2A+2C e -[A,[AC]]/3+2[C, [C,A]]/3+- _ ( B2 ) 

We write (X a )'(0) and (X a )"(0) the two first time derivatives at t = of the free-evolved observable 
X a , see ©, and similarly for B. Let us take A = -j3H v /2 and C = -(3(H B + H VB )/2. Then 

[A, [AC]] = ^-(XT(O)B (B3) 

[C, [C, A]] = -*£■ ((r)'(0) B'(0) - ^^ 2Q " 2 ^ 2 ) • (B4) 

Each time derivative of X a (of -B) gives a extra factor of T^, (tg ). Therefore, the right-hand 
side of (|B3|) is smaller than (3Hp B = PX a B by a factor of the order of (h(3/Tp) 2 <C 1. By virtue 
of A th = {j3V" (<d))- 1 / 2 and T v = (M/V"^)) 1 / 2 , the right-hand side of (EU) is of the order of 
(H(3) 2 (Tp tg + TZ 2r q th )(3Hpg with r/th given by ((ID]) . Assuming that 77th is at most of order 1 (this 
is the case in particular if (]40p holds true) this indicates that the double commutators (JB3IIB4P are 
much smaller than C and can be neglected in (JB2I) when h/3 <C T-p and (h(5) 2 <C T-p tg. Neglecting 
these commutators, (1B2|) reduces to (117p for the aforementioned choices of A and C. Notice that 
our approximation of p^H is self-adjoint and is better than Z^ >l3 e~ l3Hv e~ l3 ( HB+H ' PB ' (the error is 
of one order smaller in Hf3/Tp). The approximation e~ ltHv ' p^e v ' ~ ppg for t -C T-p is 
obtained similarly, by using (jBip with j4 = —itHp/h and C = —(3(Hp + PTg + Hpg). One can 
check explicitly by means of similar arguments as in Sec. IVII CI that the relative errors are small. 
Indeed, one can show that if one multiplies (JB3P or (jB4p by the approximate equilibrium state 



^ra e_ ' ?/fp e ~ ^ B+fip8 ' e_ ' 3Jj ' ' traces out the bath variables and takes the matrix elements 
between \x) and \x'), the matrix elements so obtained are much smaller than (]79p if ()40p is satisfied 
and h/3 <^Tp. 
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APPENDIX C: PROPERTIES OF THE BATH CORRELATION FUNCTION 

In this appendix, we establish some general properties of the bath two-point autocorrelation 
function 

h(h,t 2 ) = (B(h)B(t 2 )) = h(h - t 2 ) (CI) 

and its Fourier transform 

/oo 
dth(t)e iujt . (C2) 

-oo 

Most (but perhaps not all) of these properties are well known. The average (•) in (|Cip is taken 
with respect to the Gibbs state p B , B(t) is the bath coupling agent in the interaction picture, 
and (B) = 0, see (plfTTj), The fact that h{ti,t 2 ) depends only on the time difference t\ — t 2 is a 
consequence of the stationarity of pg [501 ] . 

Real and imaginary parts. The real and imaginary parts of h(t) are given by ?R.h(t) = {B{t)B + 
BB(t))/2 and Qh(t) = -i([B(t),B]}/2. We write (Wi)(lo) and (3fr)(w) their Fourier transforms. 
Then 

sR/t(t) = m(-t) , %h(t) = -Qh(-t) 

— ,—. ^ . — (C3) 

(m)ka) = (9th)(-u) , ($Sh)(w) = -(3f/i)(-w) . 

The imaginary part Qh(t) is linked to the linear susceptibility by x(t) = —26(t)Qh(t)/h, where 9{t) 
denotes the Heaviside function [50(]. Such a susceptibility characterizes the response of the bath 
when its Hamiltonian is perturbed by the time-dependent potential Vs(t) = —x a (t)B, where x a {t) 
is a real- valued function of time. More precisely, if B(t) is the observable B in the Heisenberg 
picture (i.e., dB(t)/dt = (i/H)[H B + V B (t),B(t)]) then (B(t)) = JdTx(t)x a (t - r) up to terms of 
order x 2a . 

The function h(t) is of positive type. This means that for any integer n > 1, complex numbers 
c\ , . . . , c n , and times t\, . . . , t n , one has 

n 

Y, c* Cj h(ti - tj) > . (C4) 

This property can be easily checked on (JC1|) . It is equivalent to h(u) > for any real u. The real 
part of h(t) is also of positive type, as 2(JR.h)(uj) = h(oS) + h(—uj) > for any real u. By using the 
Cauchy-Schwarz inequality for the Hermitian sesquilinear form (A,B) i— ► (A'B), one shows that 
\h(t)\ < h(0) = (B 2 ) for any time t. 

44 



KMS property. This property is specific to our choice of the Gibbs state for the bath average. 
It says that h(t) can be extended to an analytic function in the strip {z £ C; —H/3 < z < 0}, 



continuous on {z £ C; — hf3 < z < 0}, and such that [47] 



h(t) = h(-t - \h(5) , t € R . (C5) 

Deforming the path of integration in (|C2p . one can show that ()C5p is equivalent to h(u) = 
e^^h^—uj). In view of (|C3|) . this means that 

TOM = 'WPi* • (C6) 

By replacing in this equation (5ft/t)(o») and ($sh)(uj) by their Fourier integrals, expanding and 
identifying each power of w, and using the parity properties (|C3[) . one finds relations between the 
integrals J dtt a h(t) for even and odd a's. For instance, the identification of the zero-th power in 
uj in (|C6|) yields 

dt t h(t) = -M r dt h(t) . (C7) 

2 J — oo 

We now assume that \(^sh){uj) ~ 7Lj m as uj — > with m. a positive odd integer and 7 > (such 
a choice is motivated in Sec. IVICp . By ()C6p . this entails ($lh)(u>) ~ 27a;" 1-1 (a/3) -1 . Let a be a 
nonnegative integer, a < m — 2. Then 

dtf *W) = lim f ^(i)( W ) fdtfe-^ = -(-i^a! /°° ^^^ . (C8) 

e ^ +y-oo 27 r Jo J-00 27r ^ + 

Note that the frequency integral converges for a < m — 2, vanishes for even a's, a < m — 2, and 
diverges for a > m — 1. For a = m = 1, a similar formula holds, 

ro ° dco (Wi)(u) - 0i)(O) 

J-00 



dttm(t) = - —± — ;v y / — ^ (C9) 

27T UJ Z 



where the diverging frequency integral has been regularized by subtracting (JR.h)(0) to (8lh)(uj). 
This is equivalent to subtracting (Jfth)(0) 5(t) from $th(t) in (|C8P and this does affect the left-hand 
side of this equation. The integral on the right-hand side of (JC9P converges since (Jfth)(ui) — (Jfth)(0) 
behaves like u> 2 as u — ► 0. 

Integration in the complex plane. Let us denote by j(t) = j(—t) the primitive of Qh(t) vanishing 
at t — ► ±00, 



/ /-00 /"00 J 

dr9 ! /i(r) = — I dr 9/1(7") = / — 7(0;) cos (wi) . 
-00 Jt jo ^ 



(CIO) 
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The KMS property (|C6P and the bound tanh(u) < u for u > imply 

»sW-^s>H. (oil) 

Substituting (jCllj) into (IClOj) we obtain the following inequalities for 70 = 7(0) 

< 70 < ^h(0) . (C12) 

We can exploit the KMS property further to obtain the two identities 

»h/3 



dzh(-iz-t) = 27(f) , (C13) 

/o 

l-h/3 rz\ rh/3 

/ dzi dz 2 h{-\z 2 ) = / dzh{-\z)z = hp-fQ- ( C14 ) 

Jo Jo Jo 

To show (|C13p , we deform the contour of integration in the complex plane to get 

/•—ih/3 i-oo 

/ dz h(z - t) = / dr (h(r - t) - h(r - t - ih/3)) (C15) 

Jo Jo 

and then use (|05[) . (|C3|) . and (|C10p . The second equality in (|CJ14|) is established in a similar way, 
relying also on (|C7|) . Finally, we note that the left-hand side of (|C14[) reads 

h(3 / dz 2 /i(-iz 2 ) - / dz 2 h(-iz 2 ) z 2 . (C16) 

Jo Jo 

Therefore, the first equality in (|C14[) is a consequence of the second one and of (|C13|) . 



APPENDIX D: WICK THEOREM 



We show in this appendix that Wick's theorem (|12j) implies formulas (|70p and (|82p . i.e., 
*i,y[M = (^t,o[fe] t ^,o[/]) y exp|-iy dr(A;(r)-/(r))(5(r)),| 

= expj-^ / dn r 1 dr 2 (A;(ri) - /(n)) (fcfaj/ifa.ri) - Z(^)Mri,r 2 ))} (Dl) 

where t and y are real numbers (time and position), k and / are (nice) real functions, (-) y is the 
bath average ([83"p . h(ri,T 2 ) is the bath function (jClh . and 

C/ t ,oW = Texp{-^ I drk(T)B{r)\ . (D2) 

Let us first recall that Wick's theorem (|12p can be rephrased as the following recursive relation 
for the bath n-point functions (|10p . n > 2, 



h n {ti, . . . ,t n ) = ^2 h(t j: t i )h n ^ 2 (ti,... 1 t n )+ ^2 h(ti,tj)hn-2(ti,---,t n ) (D3) 
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wherein i is a fixed integer between 1 and n and the time arguments of the (n — 2)-point functions 
are the times appearing in the n-point function on the left-hand side except for t{ and tj. Wick's 
theorem (|D3[) holds for any value of N if the bath consists of N harmonic oscillators linearly 
coupled to the pointer (Sec. IVIDJ) . For the more general baths considered in this work, its validity 
relies on the limit N S> 1 and is a consequence of the additivity of the bath coupling agent B 
and Hamiltonian Hq in single-degree-of- freedom contributions and of the QCLT of Ref. 22|. This 



theorem provides a mapping between the correlation functions h n (ti, . . . ,t n ) and the correlation 
functions of a certain bath of harmonic oscillators in thermal equilibrium. In such a mapping, B v 



is identified with the position of the z^th oscillator [22] . 

We now proceed to proving ()Dip . Let us first consider the case y = 0. To shorten the notation, 
we write dk(r) in place of /c(r)dT. The two members of (|D1|) being equal at time t = 0, it is 
enough to prove that they satisfy the same first-order time differential equation. Hence, we need 
to show that F t: o[k,l] = (Utfllk]^ ^t,o[^]} satisfies 

dF t , [k,l\ 1 



dt h 2 



(k(t) - l(t)) / (dfc(r) h(r, t) - d/(r) h(t, T))F tfi [k, I] . (D4) 







But ihdU t fl[k]/dt = k(t)B(t)U t)0 [k], hence flD4]) is equivalent to 

ft 
U tfl [k]^B(t)U tfi [l}) =± / (dk(T)h(T,t)-dl(T)h(t,T))(u tfi [k?U t ,o[l}) . (D5) 



To show (JD5J), let us expand the two time-ordered exponentials on the left-hand side of (JD5J). 
Invoking also (|D3[) . (B(t))o = 0, and setting ho = 1, this left-hand side reads 

N 

1 

ftn+r. 



\ n (—\) m I' I' 

E ~T^r / dfc ( T l) • • • dk ( T n) / dl (h) ■ ■ ■ <U(tr. 

n+m >l J0<Tn<-<n<t J0<t m <-<ti<t 



(1 — <^n,o) 2_j h( T p>t)hn+m-l{T~n, ■ ■ ■ ,T p +l,T p -i, . . . , T\, t\, . . . ,t m ) + 
p=l 

m \ 

(1 - Sm^^j^^^n+m-liTn, . . . ,Tl,tl,... ,T q -i,T q+ \, . . . , t m ) J . (D6) 

9=1 / 

One may perform variable substitutions in the integrals in such a way that the time arguments of 
the first {n + m — l)-point function become (r n _i , . . . , n, ii, . . . t m ) and those of the second become 
(r n , . . . ,n, ti, . . . , t m -\). Doing so and resuming the series, we find that (|D6|) reduces to (|D5p . 
Turning to the case y ^ 0, we invoke (|74p to write 

U tfi [k\ ] U t>0 [l}) = ^ (U-m^m UtM ] U t ,o[l]) ( D7 ) 

I y Zj y \ / 
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wherein U-m/3fi[y a } is obtained by choosing the complex time —ih/3 and the constant function y a in 
(|D2h . By expanding the three time-ordered exponentials in the right-hand side of (|D7|) and using 
Wick's theorem (|D3|h one obtains in a similar way as above 



^ t (Ut,o[^U t ,o[l]) y = -^{k(t)-l(t))(j\dk(r)h(r,t)-dl(r)h(t,r)) + 

+\y a f dzh{-iz,A(u tfi {k^Ut, Q {l]) y ■ (D8) 



Therefore, the functional 

Gt,y[k,l] = (y tfi [k] Ji Ut,o[l\) ex P<[|r / dr / dz(fc(r)-/(r))/ J (-iz,r)| (D9) 

satisfies the same time differential equation (|D4p as .F^offc, /]. Moreover, it is equal to 1 for t = 0. 
Thus Gt i?y [fc, I] = i^o[^M] f° r an y *) 2/> fci an d ^ Setting A; = in this equation and differentiating 
with respect to 1(t) at Z = yields 

- *£ £ dzh(-iz,r) - j:(B(t)) v = ~(B(r)) = . (D10) 

This identity and (|C13() imply ([83]) . We may now replace (|D10|) into (|D9|) to get (|D1[) . Let us stress 
that, although (|D1|) and (|D10|) coincide with the lowest-order results of perturbative expansions in 
k and I, these formulas are in fact valid to all orders in k and I in the limit N ^> 1, as a consequence 
of the QCLT. 
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